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ABSTRACT 



i We have studied the optical spectra of a sample of 28 O- and early B-type stars in the Large Magellanic Cloud, 22 of which are associated with 
f~ ^ . the young star forming region Nil. Our observations sample the central associations of LH9 and LH10, and the surrounding regions. Stellar 
' parameters are determined using an automated fitting method (Mokiem et al. 2005), which combines the stellar atmosphere code fastwind (Puis 
. et al. 2005) with the genetic algorithm based optimisation routine pikaia (Charbonneau 1995). We derive an age of 7.0 ±1.0 and 3.0 ± 1.0 Myr 
• • for LH9 and LH10, respectively. The age difference and relative distance of the associations are consistent with a sequential star formation 
. scenar io in which stellar activity in LH9 triggered the formation of LH10. Our sample contains four stars of spectral type 02. From helium and 
, hydrogen line fitting we find the hottest three of these stars to be ~49 - 54 kK (compared to ~45 - 46 kK for 03 stars) . Detailed determination of 
» I ■ the helium mass fraction reveals that the masses of helium enriched dwarfs and giants derived in our spectroscopic analysis are systematically 
' lower than those implied by non-rotating evolutionary tracks. We interpret this as evidence for efficient rotationally enhanced mixing leading to 
the surfacing of primary helium and to an increase of the stellar luminosity. This result is consistent with findings for SMC stars by Mokiem et 
al. (2006). For bright giants and supergiants no such mass discrepancy is found; these stars therefore appear to follow tracks of modestly or non- 
rotating objects. The set of programme stars was sufficiently large to establish the mass loss rates of OB stars in this Z ~ 1/2 Zq environment 
sufficiently accurate to allow for a quantitative comparison with similar objects in the Galaxy and the SMC. The mass loss properties are found 
to be intermediate to massive stars in the Galaxy and SMC. Comparing the derived modified wind momenta D mom as a function of luminosity 
with predictions for LMC metallicities by Vink et al. (2001) yields good agreement in the entire luminosity range that was investigated, i.e. 
5.0 < logL/LQ < 6.1. 

Key words. Magellanic Clouds - stars:atmospheres - stars: early-type - stars: fundamental parameters - stars: mass loss 



1. Introduction 

Massive stars play an intricate role in the evolution of galax- 
ies. Because of the large energies associated with their stel- 
lar winds, ionising radiation, and life-ending supernova ex- 
plosions, they dictate galactic structuring processes such as 
star for mation and the creation and evolution of supperbub- 
bles (e.g. Oevlll999h . Mounting evidence also points to a direct 
link betwe en massive stars an d exotic phenomena such as y-ray 
bursts (e.g. Hiorth et al]|2003 ) and the reionisation of the early 



universe dBromm et al.ll200ll) . Accordingly, understanding the 
properties of these stars, both in terms of their fundamental pa- 
rameters as well as their evolution, is fundamental. 

The initial metal composition (Z) of the gas out of which 
massive stars form has a strong impact on their global prop- 
erties and characteristics. Many studies have shown that pa- 
rameters such as the effective temperature and ionising fluxes 



are strong functions of Z (e.g. Kudritzki 2002; Mokie m et al 



2004; M assev et al . 2005; M okiem et al. 2006), adding an ex- 



tra dimension to the conversion of morphological properties 
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such as spectral type to phy sical quantities. Theoret i cal and ob- 
serva tional arguments (e.g. lKudritzki & Puis Vink et al 



2001) also point to a relation between the strengths of the stel- 
lar winds of these objects and metallicity. As wind mass loss 
leads to partial evaporation of the star (e.g. Chiosi & Maeder 
19861) - with possible consequences for the nature of the com- 
pact object that is left behind after the final s upernova explosion 
- and to loss of angular momentum (e.g. IMeynet & Maedei 
2000), it also dictates to a large extent its evolutionary path 



and fate (e.g. lYoon & Langerl 120051: IWooslev & Hegerl l2006). 
Quantifying the mass loss versus metallicity dependence M(Z), 
therefore, is an important quest in astrophysics. 

Due to their proximity and low metal content the 
Magellanic Clouds provide us with unparallelled laboratories 
to test and enlarge our knowledge of massive stars. These 
galaxies, therefore, have been in the focal point of many stud- 
ies analysing the i r massive star conten t . Early studies (e.g . 



Conti etal.1 1 19861: iGarmanv etal 
Parker et al 



1987t iMassey et all 1 19891: 



19991) predominantly relied 



19921: IWalborn et al 
on photometric data and spectral type calibrations. Only rel- 
atively recent the advent of large telescopes and the develop- 
ment of sophisticated stellar atmosphere models h as allowed 
for more detailed analyses of individual stars (e .g. Puis et al 



1996 : lffillier & Milled 1 999[ ICrowther et al.l2002tlBouret et al 
20031: iMartins et al.ll2004 . Though all these studies have 



tributed enormously to our understanding of massive stars, the 
samples analysed so far have been rather limited in size (a few 
objects at a time) and have been focused predominantly on ob- 
jects in the Small Magellanic Cloud (SMC) - because its metal 
deficiency is more ext reme than that of the Lar ge Magellanic 
Cloud (LMC) (but see iMassev et alj|2004l2005b . 

The limited sizes of the samples that have so far been stud- 
ied are at the root cause of the perhaps somewhat disheartening 
conclusion that - in spite of all the progress that has been made 
- we still cannot provide robust and sound answers to the ques- 
tion: what is the role of metal content, stellar winds and rota- 
tion in the evolution of massive stars? To help attack this prob- 
lem, our research group has conducted a VLT-FLAMES Survey 
of Massive Stars (see Evans et al.l 120051) . In this ESO Large 
Program the Fibre Large Array Multi-Element Spectrograph 
at the Very Large Telescope was used to obtain optical spec- 
tra of more than 50 O- and early B-type stars in the Magellanic 
Clouds. 

Here we present the homogeneous analysis, employing au- 
tomated spectral fitting methods, of a sample of 28 O-type and 
early B-type stars in the Large Magellanic Cloud; 22 targets 
from the FLAMES survey and 6 from other sources. This is 
so far the largest sample of massive stars studied in the LMC 
and it almost doubles the amount of massive objects in this 
galaxy for which parameters have been derived from quantita- 
tive spectroscopy. Specifically, we will try to establish the mass 
loss rates of OB stars in this Z ~ 1/2 Zq environment to a level 
of precision that allows for a quantitative comparison with sim- 
ilar objects in the SMC and our Galaxy. This will provide a new 
Z-point in testing the fundamental prediction provided by radi- 
ation driven wind theory for t he mass-loss - metalli c ity depen- 
dence: M(Z) oc z - 5 "" - 7 (e.g. iKudritzki et al.lll987t IPuls et al.l 
2000; IVink etZhoOll) . 



The majority of our LMC sample is associated with 
the spectacular star forming region Nil (iHenizd Il956l) . It 
has a Ha luminosity only surpassed by that of 30 Doradus 
(Ken nicutt & Hodgd 1 19861) . ranking it as the second largest 
Hn region in the Magellanic Clouds. Nil is host to several 
OB associations of apparently different ages, the formation of 
which is believed to ha ve been triggered by stellar activity in 
the central OB cluster dParkeretal.il 19921: IWalborn & Parker! 
19921; IWalbornet ail 1999). Our observations sample both the 



central cluster LH9 as well as the younger cluster LH10, allow- 
ing for an investigation of a possible sequential star formation 
scenario. We will use our Nil stars to test predictions of mas- 
sive star evolution, including the role of rotation, and the star 
formation history. 

This paper is organised as follows: in Sect. |2]we describe 
the LMC data set that was analysed using our automated ge- 
netic algorithm based fitting method. A short description of this 
method is given in Sect.|3]and the results obtained are presented 
in Sect. H] with fits and comments on individual objects given 
in the appendix. In Sect. [5] we investigate the discrepancy be- 
tween spectroscopically determined masses and those derived 
from evolutionary tracks. The evolutionary status of Nl 1 is dis- 
cussed in Sect. [6] Finally, Sect.[7]summarises and lists our most 
important findings. 



2. Data description 

Our OB-type star sample is mainly drawn from the targets ob- 
served in the LMC within the context of the VLT-FLAMES 



survey of massive stars (see lEvans et al.1 12005). Two fields in 
the LMC centred on the clusters Nil and NGC 2004 were 
observed in the survey. Here we analyse a subset of the ob- 
jects observed in the Nil field. This set consists of all O-type 
spectra obtained, excluding those that correspond to confirmed 
binaries, and five early B-type spectra of luminous giant and 
supergiant stars. 

To improve the sampling in luminosity and temperature, we 
supplemented the FLAMES targets with six r elatively bright O- 
type field stars. These objects are p art of the San duleak (Il970i 



hereafter Sk) and lBrunet et al.l dl975l hereafter BI) catalogues, 



and were observed as part of the programs 67.D-0238, 70. D- 
0164 and 074.D-0109 (P.I. Crowther) using the Ultraviolet and 
Visual Echelle Spectrograph (UVES) at the VLT. 

The obser vations of the FLA MES targets are described ex- 
tensively by lEvans et al. (2006), to which we refer for full de- 
tails. Here we only summarise the most important observa- 
tional parameters. Basic observational properties of the pro- 
gramme stars together with common aliases are given in Tab.Q] 
Note that Nl 1-031, BI 237, BI 253 and Sk-67 166 were stud- 
ied recently using line blanketed stellar atmosphere models. In 
the appendix a comparison with these analyses is provided. The 
FLAMES targets were observed with the Giraffe spectrograph 
mounted at UT2. For six wavelength settings a spectrum was 
acquired six times for each object with an effective resolving 
power of R ^ 20 000. These multiple exposures, often at differ- 
ent epochs, allowed for the detection of variable radial veloc- 
ities. As a result, a considerable number of binaries could be 



M. R. Mokiem et al.: Wind properties and evolution of hot massive stars in the LMC 



3 



Table 1. Basic parameters. Primary identificatio n numbers for Ni l are from Evans et al. I (120061). Identifications starting with 
"Sk", "PGMW" and "BI" are from, respectively ISanduleakl dl970l) IParker et all d 19921) and iBrunet etall d 1975b . Photometric 
data for these objects are from|Evans et al. d2006b and IParker et alJ d 19921). t he latt er are flagged with an asterisk. For non-Nl 1 
object s these are fromlArdeberg et a" ( 1972irisserstedl (1 1979b andlMassevI d2002b . Wind velocities given without brackets are 
from lCrowther et alj d2002h . iMassa et al] d2003b and lMassev et alj d2005b . For Sk -66 18 the wind velocity was measured from 
Oiv (1031-1037 A) in its FUSE spec trum. Values between brackets are calculated from the escape velocity at the stellar surface. 
F or Nl 1-03 1 the value of v x is from lWalborn et al] (12004b . though this is actually based on the value obtained for Sk-68 137 by 



Prinja & Crowt her ( 1998). 



Primary ID 


Cross-IDs 


Spectral 
Type 


V 


A v 


M v 


Woo 

[km s~'] 


Nl 1-004 


Sk -66 16 


OC9.7 lb 


12.56 


0.74 


-6.68 


[2387] 


Nl 1-008 


Sk -66 15 


B0.7 la 


12.77 


0.84 


-6.57 


[1619] 


Nl 1-026 




02 III(f*) 


13.51 


0.47 


-5.46 


[3116] 


Nl 1-029 




09.7 lb 


13.63 


0.56 


-5.43 


[1576] 


Nll-031 


PGMW 3061 


ON2 III(f*) 


13.68* 


0.96 


-5.78 


3200 


Nl 1-032 


PGMW 3168 


07 11(f) 


13.68* 


0.65 


-5.47 


[1917] 


Nl 1-033 


PGMW 1005 


B0 Illn 


13.68 


0.43 


-5.25 


[1536] 


Nl 1-036 




B0.5 lb 


13.72 


0.40 


-5.18 


[1714] 


Nl 1-038 


PGMW 3100 


05 n(f+) 


13.81* 


0.99 


-5.68 


[2601] 


Nl 1-042 


PGMW 1017 


B0 III 


13.93 


0.22 


-4.79 


[2307] 


Nl 1-045 




09 III 


13.97 


0.50 


-5.03 


[1548] 


Nl 1-048 


PGMW 3204 


06.5 V((f» 


14.02* 


0.47 


-4.95 


[3790] 


Nl 1-051 




05 Vn((f)) 


14.03 


0.19 


-4.66 


[2108] 


Nl 1-058 




05.5 V((f)) 


14.16 


0.28 


-4.62 


[2472] 


Nl 1-060 


PGMW 3058 


03 V((f*)) 


14.24* 


0.81 


-5.07 


[2738] 


Nl 1-061 




09 V 


14.24 


0.78 


-5.04 


[1898] 


Nl 1-065 


PGMW 1027 


06.5 V((f)) 


14.40 


0.25 


-4.35 


[2319] 


Nl 1-066 




07 V((f» 


14.40 


0.25 


-4.35 


[2315] 


Nl 1-068 




07 V((f») 


14.55 


0.28 


-4.23 


[3030] 


Nl 1-072 




B0.2 III 


14.61 


0.09 


-3.98 


[2098] 


Nl 1-087 


PGMW 3042 


09.5 Vn 


14.76* 


0.62 


-4.36 


[3025] 


Nll-123 




09.5 V 


15.29 


0.16 


-3.37 


[2890] 


BI 237 




02 V((f*)) 


13.89 


0.62 


-5.23 


3400 


BI 253 




02 V((f*)) 


13.76 


0.71 


-5.45 


3180 


Sk -66 18 




06 V((f)) 


13.50 


0.37 


-5.37 


2200 


Sk -66 100 




06 11(f) 


13.26 


0.34 


-5.58 


2075 


Sk -67 166 


HD 269698 


04 Iaf+ 


12.27 


0.31 


-6.54 


1750 


Sk -70 69 




05 V 


13.95 


0.28 


-4.83 


2750 



detected (Eva ns et al. 2006), which we subsequently excluded 
from our analysis. 

To allow for a sky subtraction a master sky-spectrum was 
created from combining the sky fibres in the Giraffe spectro- 
graph (typically 15), individually scaled by their relative fibre 
throughput. Even though in general the sky background is low 
and this approach successfully removes the background con- 
tribution, in crowded regions such as Nl 1 accurate subtraction 
of nebular features remains very difficult. As a result of this, 



the line profiles of many of our programme stars still suffer 
from nebular contamination. This in principle does not hamper 
our analysis. For most stars the core nebular emission is well- 
resolved and we simply disregard this contaminated p art of the 
profile in the automated line fits. lMokiem et al] (2006) showed 
by performing tests using synthetic data that with this reduced 
amount of information, the automated method can still accu- 
rately determine the correct fit parameters. Also tests assessing 
the impact of possible residual nebular contamination in the 
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line wings or over-subtraction of sky components showed that 
its effect is negligible. 

For each wavelength range the individual sky-subtracted 
spectra were co-added and then normalised using a cubic- 
spline fit to the continuum. A final spectrum covering 3850- 
4750 and 6300-6700 A was obtained by merging the nor- 
malised data. Depending on the magnitude of the target these 
combined spectra have typical signal-to-noise ratios of 100- 
400. 

The first four field stars listed in Tab. Q]were observed with 
the VLT-UVES spectrograph in service mode on 29 and 30 
November 2004 under program 74.D-0109. UVES is a two 
armed cross-dispersed echelle spectrograph allowing for simul- 
taneous observations in the blue and red part of the spectrum. 
In the blue arm standard settings with central wavelengths of 
390 and 437 nm were used to observe the spectral ranges 3300- 
4500 and 3730-5000 A. For the red arm standard settings with 
central wavelengths of 564 and 860 nm provided coverage be- 
tween 4620-5600 and 6600-10400 A. A 1.2" wide slit was 
used, providing a spectral resolution of 0.1 A at Hy, corre- 
sponding to an effective resolving power of R ^ 40 000, a value 
which applies to all UVES setups. Individual exposure times 
ranged from 1500 to 2200 seconds. 

Sk -67 166 was observed on 27 and 29 September 2001 
under program 67.D-0238 with UVES using a 1" slit. A stan- 
dard blue setting with central wavelength 437 nm provided con- 
tinuous coverage between 3730-5000 A. A non-standard red 
setup with central wavelength 830 nm was used to observe 
the range between 6370-10250 A. Three exposures, each of 
1000 s were obtained for each setup. Note that the data fo r 
this target was previously presented bv lCrowther et al. (2002). 
Finally, Sk -70 69 was observed on 1 and 2 December 2002 
using UVES under program 70.D-0164. Standard blue and red 
settings with central wavelengths 390 and 564 nm were used 
in simultaneous 2400 s exposures using dichroic. A second 
non-standard red setup with central wavelength 520 nm (4170- 
6210 A) without dichroic was used in a 1500 s exposure. The 
typical two pixel S/N ratios obtained for all spectra are 80 at 
Hy and 60 at Ha. 

Spectral types for FLAMES stars were determined by vi- 
sual inspection of th e spectra, using published sta ndards. In 
particular the atlas of IWalborn & Fitzpatrick (1990) was used 
while considering the lower metallicity environ ment of the 



LMC. The classifications given in Tab.Q]are from Evans et al 
d2006l) . A comparison with previously published spectral types 



is also given in the latter pap er. For the field s t ars we adopted 



the classifica tions given by IWalbornl d!977l) . 
dl995il2002al) and lMassev et all d 19951 l2005l) . 



Walborn et al 



Photometric data for the FLAMES targets was obtained 
predominantly from B and V images of the Nil field 
taken with the Wide Field Imager at the 2.2-m Max Planck 



Gesel lschaft/ESO telescope on 2003 April 2004 dEvans et al 
2006). The p hotometry for stars flagged with an asterisk was 



adopted from [Parker et al.1 (119921) . The caption of Tab. Q] lists 
the references to the various sources of the photometric data of 
the field stars. 

To calculate the interstellar ext inction (Ay) giv en in Tab.Q] 
we adopted intrinsic colours from Johnson ( 19661 and refer- 



ences therein) and a ratio of total to selective extinction of 
Rv = 3.1. With these Ay values, extinction corrected visual 
magnitudes (Vb) were calculated from the observed V-band 
magnitudes. Finally, we calculated the absolute visual magni- 
tude My, while adopting a distance modulus of 18.5 for the 



tude My, while adopting a distance modulus o 
LMC dPanagia et al!l99lHMitchell et al.ll2002h . 



3. Analysis method 

All optical spectra ar e analysed using the automated fitting 
meth od devel oped by iMokiem et al.l ((2005, hereafter referred 
to as IPaper I ). Here we will suffice with a short description 
of the method and re fer to the before m entioned paper and to 
Mokiem et al.l (120061 hereafter IPaper III) for the details. 

In short the automated fitting metho d uses the genetic algo- 
rithm based optimisation routine pikaia dCharbonneaul 1 9951) to 
determine the s et of input parame ters for the stellar atmosphere 
COde FASTWIND dPulsetal.ll2005l) which fit an observed spec- 
trum the best. This best fitting model is constructed by evolving 
a population of fastwind models over a course of generations. 
At the end of every generation the parameters of the models 
which relatively fit the observed spectrum the best are used to 
construct a new population of models. By repeating this proce- 
dure a natural optimisation is obtained and after a number of 
generations (see below) the best fitting model, i.e. global opti- 
mum in parameter space, is found. 

Using the concept of a unified model atmosphere the fast 
performance code fastwind incorporates non-LTE and an ap- 
proximate approach to line blanketing to synthesise hydrogen 
and helium line profiles. Consequently, given the spectral range 
of the observed data set in this study we will focus on the mod- 
elling of the optical hydrogen and helium lines. To account for 
the accuracy with which each individual line can be reproduced 
b y fastwin d we adopt the line weighting scheme as described 
inlj - 



Paper I . 



In the fitting of the spectra we allow for seven free parame- 
ters. These are the effective temperature T^, the surface gravity 
g, the helium number density defined as Yn e = 7V(He)/7V(H), 
the microturbulent velocity u tul -b, the projected rotational veloc- 
ity u r sin i, the mass loss rate M and the exponent of the beta- 
type velocity law describing the supersonic regime of the stellar 
wind. For the terminal velocity of the wind 0^, which cannot 
be accurately determined from the optical spectrum, we adopt 
values determined from the analysis of ultraviolet (UV) wind 
lines. If no UV determination of » M is available a scaling re- 
lation of Woo with the escape velocity (u eS c) defined at the stel- 
lar surface is used throughout the fitting process. For our pro 



gramme stars we adopt the ratio as determined bv lLamers et al 



(1993 of foo/fesc — 2.6. Similar as in Paper II we adopt fixed 
values for the atmospheric abundances of the background met- 
als. These were scaled with respect to mass ratios based on 



the Solar abundances of Grevesse & Sauval (1998, and ref- 



erences therein). The metallicity scaling factor was set equal 
to the mean metal deficiency of 0.5 as foun d for the LMC 
dRussell & Besselllll989l:lRolleston et al.ll2002h . 

For our current data set the spectral range and quality, with 
exception of the signal-to-noise ratio, is similar to the set anal- 
ysed in IPaper III Consequently, we adopt the same minimum 
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number of generations that have to be calculated to determine 
the best fit, i.e. t o assure t hat the global optimum in parameter 
space is found. In lPaper fl , using formal tests that accounted for 
nebular contamination of the line profiles and a signal-to-noise 
ratio of 50, this number was determined to be 150. 

The uncertainties of the fit parameters are determin ed us- 
ing so-called optimum width based error estimates. In lPaper II 
we have argued that an error estimate for a parameter can be 
defined as the maximum variation of this parameter within the 
global optimum in parameter space. By measuring the width 
of this optimum in terms of fit quality one can determine which 
models, based on their fit quality, are associated with the global 
optimum. The maximum variation of the individual fit param- 
eters within th is group of models then gives the error estimates 
(see Sect. 4 in lPaper I ). 

In Tab.[3]the optimum width based error estimates of the fit 
parameters for our programme stars are listed. The uncertain- 
ties in R+, L+, M s and M ev , that are derived from the fit param- 
eters, were calculated using the same approach as in iPaper II . 
The adopted uncerta inty in the visual magnitude is AMy = 0.13 
(Pana giaet alJll99ll) for all our programme stars. 



3.1. The assumption of spherical symmetry 

fastwind assumes a spherically symmetric star and wind. For 
very high rotational speeds, u ml , this may potentially be a prob- 
lematic assumption. As it is well known, a high rotation rate 
leads to a distortion of the stellar surface and, via the von Zeipel 
theorem, to a decrease i n flux and effectiv e temperature from 
the pole to the equator (I von Zeipellll924l (slightly) modified 
by Maeder (1999) for the relevant case of shellular, i.e., radi- 
ally dependent, differential rotation, u> = u>{r)). This so-called 
"gravity" darkening d oes not only affect the stellar parameters, 
but also the wind (e.g. ICranmer & Owockill995 ; Owocki et al 



199(§ IPetrenz & Pulsll2000h . Thus, it might be questioned in 
how far the derived properties (which then depend on incli- 
nation angle and have a somewhat local character) are repre- 
sentative for the global quantities (e.g., mass, luminosity and 
mass-loss rate) which refer to integral quantities. 

A s shown by various simulations ranme r & Owockil 

19951: IPetrenz & Pulsll996ll2000b . the difference between local 



and global quantities remains small unless the star rotates faster 
than <: 60. . .70 % of its critical speed. Unfortunately, however, 
we cannot directly access the actual rotational speed, v Iot , but 
only its projected value, u r sin i. 

Accounting for the average value of < sin i >= tt/4, we 
have calculated the ratio of average to critical rotational speed, 
Q as v T sin i/(< sin i > u CI j t ), and found that only two objects 
lie above this value, namely Nl 1-033 and Nl 1-051, both with 
Q w 0.7. All other objects lie well below Q. = 0.3 (except for 
Nl 1-087 with Q. « 0.4). Of course, we cannot exclude that also 
the latter objects lie above the decisive threshold (if observed 
pole-on), but the rather large sample size implies that such a 
possibility should be actually present only for a minority of ob- 
jects. In conclusion, we suggest that the majority of our objects 
is, if at all, only weakly affected by a distortion of surface and 
wind. Note, e.g., that Q = 0.3 leads to a difference in T e ff and 



R+ between pole and equator of less than 5%, i.e., of the same 
order or less than our error estimates (cf. Tab. [3). Regarding the 
derived gravities and masses, finally, we have applied a consis- 
tent centrifugal correction anyhow (cf. Sect. l4.3b . 

For the two objects with large rotational speeds, on the 
other hand, it is rather possible that we observe them al- 
most equator-on, i.e., the observed profiles contain an intrin- 
sic averaging over the complete stellar disk and thus corre- 
spon d, at least in part, t o the global, polar angle averaged val- 



ues. |Howarfli&_Smith (2001) analyzed two galactic fast rota- 
tors, HD 149757 (^Oph) and HD 191423, accounting for ef- 
fects of non-sphericity and von Zeipel's theorem. Their study 
showed that for these objects, rotating at Q. = 0.9, differ- 
ences in effective temperature and radius between pole and 
equator can am ount to 20-30 %. Int e restingly, analysis of th e 
sam e objects by | Herrero et aD d2002b . IVillamariz et alJ d2002h . 
and IVillamari z & Herrer ol d2005lT using both hydrostatic and 
spherically symmetric expanding atmospheres yielded (aver- 
age) parameters that agreed to within the standard devi ation 
(5-10%) with the results obtained by iHowarfh & Smith! For 
stellar parameters it thus appears justified to conclude that non- 
sphericity has a very small impact even on stars with extreme 
v t sin i. But note also that with respect to global mass-loss rates 
the situation is much more unsecure and more-D simulations 
would be required to constra in their actual values. As shown, 
e.g., by IPetrenz & Puis (2000), the modified wind-momentum 
rate from a rapidly rotating (Q = 0.85) B-supergiant might be 
underestimated up to one magnitude if seen equator on. As the 
wind properties of Nl 1-031 and Nl 1-051 are rather normal it 
seems however unlikely that in these cases rotational effects 
have this type of dramatic impact. 



4. Fundamental parameters 

4.1. Effective temperatures 

In Fig. Q] the distribution of the effective temperatures of our 
programme stars as a function of spectral type is shown. The 
different luminosity classes are denoted using circles, trian- 
gles and squa res, respe ctively, for class V, III and I-II objects. 
Similar as in IPaper III we see that for a given spectral type 
the dwarfs are systematically hotter than the giants and su- 
pergiants. This separation can be interpreted as the result of 
the reduced surface gravities of the more evolved objects. A 
low er surface gravity re sults in an increased helium ionisation 
fe.g. lKudritzki et al J 1983b. reducing th e T e ff needed for a given 
spectral type (e.g. lMokiem et al ]l2004b . A second reason is that 
the more evolved objects have stronger winds. These denser 
winds induce an increased line bla nketing effect, further reduc - 
ing the required tem perature (e.g. lSchaerer & de Koterlll997h . 
Massey et al.1 (I2005r) . who also analysed a sizeable sample of 
LMC stars, do not find evidence for the dwarfs being hotter 
than the giants. Though, we note that their analysis only con- 
tained two giant T e ff determinations for spectral type later than 
02. 

For comparison we also show in Fig. Q] as a dashed line 
the observe d r e ff vs. spec t ral cal ibration for Galactic O-type 
dwarfs from Martins et al. (2005a). With a dotted line the av- 
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Table 2. Fundamental parameters of the LMC sample determined using GA optimised spectral fits, with T e ff in kK, log g and 
\ogg c in cms -2 , R+ in Rq, L* in L©, u tur b and u r sin/ in kms -1 , M in Mgyr 'and M s and M ev in Mq. Results were obtained 
using a population of 72 fastwind models evolved over a minimum of 150 generations. Gravities corrected for centrifugal 
acceleratio n (logg c ) were use d to calculate the spectroscopic masses (M s ). Evolutionary masses (M ev ) were derived from the 
tracks of ISchaerer etail dl993l) . 



ID 


ST 




log g 


log 9c 


R* 


logL* 




^turb 


v Y sin i 


M 




P 


M s 




Nl 1-004 


OC9.7 lb 


31.6 


3.36 


3.37 


26.5 


5.80 


0.10 


5.7 


81 


1. /O 


1 A" 

1U 




1.18 


59.9 


47.5 


Nl 1-008 


B0.7 la 


26.0 


2.98 


2.99 


29.6 


5.55 


0.10 


18.0 


83 


A OA 


1U 




1.87 


31.4 


31.7 


Nl 1-026 


02 III(f*) 


53.3 


4.00 


4.00 


10.7 


5.92 


0.11 


19.0 


109 


1 G 1 
1,51 


1U 




1.08 


42.4 


81.9 


Nl 1-029 


09.7 lb 


29.4 


3.23 


3.24 


15.7 


5.21 


0.07 


19.1 


77 


1. / J 


1U 


-7 


1.63 


15.5 


24.9 


Nll-031 


ON2 III(f*) 


45.0 


3.85 


3.86 


13.7 


5.84 


0.10 


20.0 


116 




1U 


-6 


0.89 


49.3 


60.8 


Nl 1-032 


07 11(f) 


35.2 


3.45 


3.46 


14.0 


5.43 


0.09 


11.3 


96 


£ OA 
o.UO 


1 fV 

1 u 


-7 


1.03 


20.6 


33.8 


Nl 1-033 


B0 Illn 


27.2 


3.21 


3.35 


15.6 


5.07 


0.08 


17.9 


256 


AA 


1 fV 

1 u 




1.03 


19.8 


21.0 


Nl 1-036 


B0.5 lb 


26.3 


3.31 


3.32 


15.6 


5.02 


0.08 


13.6 


54 


1 HA 
1 .UO 


1 fV 

1 u 




0.80 a) 


18.4 


20.0 


Nl 1-038 


05 II(f+) 


41.0 


3.72 


3.74 


14.0 


5.69 


0.10 


9.2 


145 


1 . JZ 


1 (Y 

1 u 




0.98 


38.8 


48.3 


Nl 1-042 


B0 III 


30.2 


3.69 


3.70 


11.8 


5.01 


0.10 


4.0 


42 


1.07 


1U 




1.19 


25.1 


21.4 


Nl 1-045 


09III 


32.3 


3.32 


3.35 


12.0 


5.15 


0.07 


16.8 


105 




1U 




0.80 a) 


11.8 


24.6 


Nl 1-048 


06.5 V((f) 


40.7 


4.19 


4.20 


9.9 


5.38 


0.06 


1.0 


130 


1.0 / 


1U 


-7 


0.80 a) 


56.6 


36.6 


Nl 1-051 


05 Vn((f)) 


42.4 


3.75 


3.88 


8.4 


5.31 


0.08 


19.7 


333 


1.U1 


1U 


-6 


0.60 


19.5 


36.4 


Nl 1-058 


05.5 V((f) 


41.3 


3.89 


3.90 


8.4 


5.27 


0.10 


14.5 


85 




1U 


-7 


1.42 


20.3 


34.4 


Nl 1-060 


03 V((f*)) 


45.7 


3.92 


3.93 


9.7 


5.57 


0.12 


19.3 


106 


5.22 


10" 


-7 


1.26 


29.2 


49.4 


Nl 1-061 


09 V 


33.6 


3.51 


3.52 


11.7 


5.20 


0.09 


19.8 


87 


2.14 


10" 


-7 


1.80 


16.7 


26.6 


Nl 1-065 


06.5 V((f) 


41.7 


3.89 


3.90 


7.4 


5.17 


0.17 


9.4 


83 


3.63 


10" 


-7 


0.80°' 


15.8 


32.9 


Nl 1-066 


07 V((f)) 


39.3 


3.87 


3.88 


7.7 


5.10 


0.11 


4.8 


71 


4.08 


10" 


-7 


0.80°' 


16.2 


29.5 


Nl 1-068 


07 V((f)) 


39.9 


4.13 


4.13 


7.1 


5.06 


0.10 


15.6 


54 


3.43 


10" 


-7 


1.12 


25.2 


29.2 


Nl 1-072 


B0.2 III 


30.8 


3.78 


3.78 


7.9 


4.70 


0.12 


7.6 


14 


2.35 


10" 


-7 


0.84 


13.8 


17.7 


Nl 1-087 


09.5 Vn 


32.7 


4.04 


4.09 


8.9 


4.91 


0.10 


15.5 


276 


1.38 


10" 


-7 


0.80 a) 


35.6 


20.9 


Nll-123 


09.5 V 


34.8 


4.22 


4.23 


5.4 


4.58 


0.09 


9.0 


110 


7.62 


10" 


-8 


0.80 a) 


17.8 


18.8 


BI 237 


02 V((f*)) 


53.2 


4.11 


4.11 


9.7 


5.83 


0.10 


12.8 


126 


7.81 


10" 


-7 


1.26 


44.6 


75.0 


BI 253 


02 V((f*)) 


53.8 


4.18 


4.19 


10.7 


5.93 


0.09 


18.6 


191 


1.92 


10" 


-6 


1.21 


64.6 


84.1 


Sk -66 18 


06 V((f)) 


40.2 


3.76 


3.76 


12.2 


5.55 


0.14 


10.8 


82 


1.07 


10" 


-6 


0.94 


31.5 


40.7 


Sk -66 100 


06 11(f) 


39.0 


3.70 


3.71 


13.6 


5.58 


0.19 


8.7 


84 


8.81 


10" 


-7 


1.27 


34.7 


41.4 


Sk -67 166 


04 Iaf + 


40.3 


3.65 


3.66 


21.3 


6.03 


0.28 


20.0 


97 


9.28 


10" 


-6 


0.94 


75.0 


70.4 


Sk -70 69 


05 V 


43.2 


3.87 


3.88 


9.0 


5.41 


0.17 


16.1 


131 


1.03 


10" 


-6 


0.78 


22.7 


39.7 



n> assumed fixed value 



erage temperature of the SMC dwarfs studied in Paper III is 
shown. The LMC dwarfs are found to occupy the region in 
between these two average temperature scales, with an average 
behaviour intermediate to that of the Galactic and SMC dwarfs. 
We interpret this as the result of the metallicity of the LMC that 
is lower than the Galactic value and higher than in the SMC. As 
the amount of line blanketing in a stellar atmosphere depends 
on metallicity the LMC objects have temperatures in between 
that of objects in the other two galaxies. 



4.2. The r e ff scale of 02 stars 

Our sample co ntains four Q2 type st ars. This spectral type was 
introduced by lWalborn et alJ ( 2002bl) and is assigned based pri- 
marily on the ratios of se lective emission lines of N iv and N m. 
By modelling these lines lWalborn et al.l (120041) have shown that 
indeed, as had been hypothesised, the 02 spectra correspond to 
higher effective temperatures. However, the correct treatment 
of nitrogen lines in stellar atmosphere models is notoriously 
difficult. The relevant ionisation stages of this atom represent 
much more complex ion models compared to the relatively 
simple hydrogen and helium ions. Moreover, for the higher ions 
of nitrogen the ionisation depends on the extreme-UV radiation 
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Table 3. Optimum width based error estimates for the seven fit parameters. The ND entries correspond to error in v tur \, that reach 
up to the maximum allowed value of i> tul -b and, therefore, are formally not defined. See text for details on the calculation of the 
uncertainties in the derived parameters. Units: r e +r in kK, \ogg and \ogg c in cms -2 , R+ in Rq, L+ in Lq, u tur \, and v, sin / in 
kms -1 , M in M Q yr~'and M s and M ev in M©. 



ID 


Ar eff 


Alog g c 


Afl* 


A log U 




M urb 


Au r sin ( 


log AM 




AM, 


AM ev 


N 11-004 


-0.5 
+0.5 


-0.06 
+0.06 


+ 1.7 


±0.06 


-0.01 
+0.02 


-2.5 
+8.2 


-6 
+7 


-0.14 
+0.05 


-0.00 
+0.19 


-10 

+ 11 


-4 
+4 


N 11-008 


-0.6 
+0.5 


-0.06 
+0.05 


+ 1.9 


±0.07 


-0.01 
+0.04 


-4.2 
+ND 


-8 
+9 


-0.20 
+0.16 


-0.31 
+ 1.08 


-5 
+8 


-2 
+2 


Nl 1-026 


-3.9 
+0.S 


-0.05 
+0.05 


±0.8 


±0.14 


-0.01 
+0.02 


-4.8 
+ND 


-11 
+13 


-0.05 
+0.10 


-0.12 
+0.03 


-7 
+9 


-16 
+21 


N 11-029 


-0.6 
+0.S 


-0.05 
+0.06 


±1.0 


±0.07 


-0.01 
+0.02 


-3.9 
+ND 


-6 
+7 


-0.13 
+0.16 


-0.42 
+0.31 


-2 
+4 


-2 
+2 


Nll-031 


-1.6 
+2.2 


-0.06 
+0.07 


±0.9 


±0.10 


-0.02 
+0.02 


-7.6 
+ND 


-22 
+21 


-0.05 
+0.06 


-0.04 
+0.06 


-9 
+ 13 


-8 
+ 11 


Nl 1-032 


-0.7 
+0.4 


-0.07 
+0.06 


±0.9 


±0.06 


-0.02 
+0.02 


-6.1 
+6.6 


-8 
+ 11 


-0.16 
+0.09 


-0.11 
+0.28 


-4 
+5 


-2 
+2 


Nl 1-033 


-0.9 
+1.0 


-0.06 
+0.09 


±1.0 


±0.08 


-0.01 
+0.03 


-5.9 
+ND 


-211 
+20 


-0.26 
+0.10 


-0.20 
+0.74 


-3 
+5 


-2 
+2 


Nl 1-036 


-0.3 
+ 1.1 


-0.09 
+0.06 


±1.0 


±0.09 


-0.01 
+0.02 


-4.2 
+ 1.4 


-4 
+6 


-1.99 
+0.30 


- 


-4 
+6 


-2 
+2 


Nl 1-038 


-1.0 
+0.8 


-0.06 
+0.07 


±0.9 


±0.07 


-0.01 
+0.03 


-9.0 
+8.3 


-19 
+18 


-0.16 
+0.10 


-0.13 
+0.22 


-6 
+ 16 


-4 
+4 


N 11-042 


-0.8 
+0.6 


-0.07 
+0.09 


±0.7 


±0.07 


-0.01 
+0.03 


-3.8 
+4.2 


-6 
+4 


-0.47 
+0.40 


-0.47 
+0.77 


-5 
+ 14 


-1 
+2 


N 11-045 


-1.2 
+0.6 


-0.09 
+0.21 


±0.8 


±0.08 


-0.02 
+0.03 


-7.6 
+ND 


-15 
+17 


-0.20 
+0.14 


- 


-2 
+5 


-2 
+3 


N 11-048 


-2.0 
+1.6 


-0.21 
+0.06 


±0.7 


±0.10 


-0.01 
+0.02 


-0.8 
+ 13.2 


-27 
+38 


-1.48 
+0.74 


- 


-22 
+30 


-4 
+5 


Nll-051 


-1.9 
+ 1.1 


-0.06 
+0.06 


±0.6 


±0.09 


-0.02 
+0.03 


-8.4 
+ND 


-38 
+23 


-0.28 
+0.11 


-0.09 
+0.52 


-3 
+5 


-4 
+5 


Nl 1-058 


-0.5 
+0.7 


-0.06 
+0.05 


±0.5 


±0.06 


-0.01 
+0.01 


-5.9 
+3.5 


-8 
+7 


-0.36 
+0.08 


-0.20 
+0.60 


-3 
+4 


-2 
+2 


N 11-060 


-1.0 
+2.3 


-0.05 
+0.09 


±0.7 


±0.10 


-0.02 
+0.03 


-13.1 
+ND 


-16 
+ 16 


-0.22 
+0.06 


-0.21 
+0.20 


-5 
+9 


-7 
+8 


Nl 1-061 


-0.6 
+0.6 


-0.09 
+0.07 


±0.7 


±0.06 


-0.01 
+0.03 


-5.4 
+ND 


-9 
+ 11 


-0.21 
+0.14 


-0.49 
+0.38 


-4 

+6 


-2 
+2 


Nl 1-065 


-0.9 
+0.7 


-0.07 
+0.15 


±0.5 


±0.06 


-0.03 
+0.04 


-9.2 
+6.3 


-8 
+8 


-1.06 
+0.14 




-3 
+3 


-2 
+2 


Nl 1-066 


-1.3 
+1.8 


-0.15 
+0.19 


±0.5 


±0.09 


-0.03 
+0.05 


-4.5 
+ND 


-15 
+22 


-0.49 
+0.22 




-5 
+ 12 


-4 
+3 


Nl 1-068 


-1.0 
+0.8 


-0.19 
+0.07 


±0.4 


±0.07 


-0.02 
+0.03 


-14.8 
+ND 


-13 
+13 


-0.79 
+0.16 


-0.24 
+0.36 


-9 
+ 14 


-2 
+2 


Nl 1-072 


-0.6 
+0.6 


-0.07 
+0.09 


±0.5 


±0.06 


-0.02 
+0.02 


-2.8 
+2.5 


-6 
+8 


-1.29 
+0.24 


-0.32 
+0.31 


-3 
+5 


-1 
+1 


Nl 1-087 


-0.6 
+1.0 


-0.09 
+0.11 


±0.6 


±0.07 


-0.03 
+0.02 


-6.7 
+.ND 


-20 
+21 


-2.82 
+0.27 




-8 
+ 14 


-1 
+2 


Nll-123 


-0.7 
+0.6 


-0.11 
+0.00 


±0.3 


±0.06 


-0.01 
+0.02 


-5.1 
+6.6 


-11 

+ 13 


-2.03 
+0.35 




-4 
+7 


-1 
+ 1 


BI237 


-3.8 
+5.2 


-0.08 
+0.15 


±0.8 


±0.18 


-0.01 
+0.03 


-12.6 
+ND 


-15 
+23 


-0.16 
+0.12 


-0.17 
+0.16 


-9 
+22 


-18 

+25 


BI 253 


-5.5 
+5.8 


-0.15 
+0.14 


±0.9 


±0.19 


-0.02 
+0.03 


-18.0 
+ND 


-21 
+18 


-0.05 
+0.10 


-0.21 
+0.07 


-21 

+39 


-23 
+31 


Sk -66 18 


-1.1 
+0.6 


-0.13 
+0.09 


±0.8 


±0.07 


-0.03 
+0.04 


-10.4 
+6.4 


-12 
+15 


-0.05 
+0.06 


-0.10 
+0.11 


-9 
+9 


-3 
+4 


Sk -66 100 


-1.0 
+0.8 


-0.09 
+0.08 


±0.9 


±0.07 


-0.03 
+0.05 


-8.1 
+8.3 


-10 
+13 


-0.12 
+0.12 


-0.19 
+0.18 


-8 
+ 12 


-3 
+4 


Sk -67 166 


-0.8 
+0.9 


-0.08 
+0.00 


±1.3 


±0.07 


-0.04 
+0.06 


-8.0 
+ND 


-23 
+15 


-0.05 
+0.02 


-0.04 
+0.09 


-15 
+35 


-7 
+7 


Sk -70 69 


-1.4 
+ 1.0 


-0.14 
+0.13 


±0.6 


±0.08 


-0.03 
+0.05 


-13.1 
+ND 


-15 

+22 


-0.11 
+0.14 


-0.25 
+0.15 


-7 
+7 


-3 
+5 



field, which may be affected by non-thermal processes, such 
as shocks. Adding to the complexity is the fact that the nitro- 
gen abundance has to be treated as a free parameter, as many 
early type stars s how evidence of atmospheric abundance en- 
hancements (e.g. ICrowfher et aL 2002 ; Bouret et al. 2003b . In 
our analysis method we solely model the hydrogen and helium 
lines and self consistently allow for abundance enhancements 
by treating the helium abundance as a free parameter. In prin- 
ciple our temperature determination should, therefore, not be 
affected by such problems. As a result of this, our analysis can 
provide an independent confirmation for the hot nature of 02 
stars. 

In Fig. [TJ we see that the objects with an 02 spectral type 
indeed correspond to the hottest stars in our sample. With ex- 



ception of the giant Nl 1-031, which based on its helium lines 
we find to be cooler (see below), they have temperatures in ex- 
cess of 50 kK and, therefore, are significantly hotter than the 
03 star at T e ff = 46 kK. The error bars, though, are consider- 
able. Compared to the average error of ~3 percent the 02 ef- 
fective temperatures have an uncertainty of 5 up to 11 percent. 
This large uncertainty can be explained by the weak or even 
absent neutral helium lines in the spectra of these objects. As a 
result of this, our fitting method has to predominately rely on 
the He n line profiles to determine the correct helium ionisation 
equilibrium. Based on a single ionisation stage the determina- 
tion of this equilibrium is more uncertain, which explains the 
larger error bars. One should also be wary for a possible degen- 
eracy effect between r e +r and Tee that can occur as a result of 
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02 03 04 05 06 07 08 09 BO B1 
spectral type 

Fig. 1. Effective temperatures as a function of spectral type for 
the investigated LMC sample. Different luminosity classes are 
denoted using circles, triangles and squares for class V, III and 
I-II objects , respectively. Shown a s a dashed line is the T e ff cal- 
ibration of iMartins et al. ( 2005a ) for Galactic O-type dwarfs. 
The dotted line corresponds to the average effective tempera- 
ture of the SMC dwarfs studied in lPaperllt The LMC dwarfs 
are found to lie in between these two average scales, with a 
typical ~2 kK offset from both the average Galactic and SMC 
relations. 



the missing neutral helium lines (see e.g. Paper J) . This, how- 
ever, is not an issue as the helium abundances derived for all 
02 stars correspond to normal values close to 0.10. 

It is also possible to question whether the fact that the He i 
lines are so weak could result in a systematically higher T e g 
determination by the automated method compared to "by eye" 
fits. In other words, would a "by eye" fit prefer a lower tem- 
perature solution? For the 02 giant Nl 1-026 this is a relevant 
question, as the weak Hei A4411 line, shown in Fig. IA.ll in the 
appendix, is not fitted perfectly. Note, however, that due to the 
relative plot scale the discrepancy is exaggerated and that the 
fit quality is good and is comparible to that of the He n lines. 
Still, to assess whether the fit of this particular line could be 
improved, i.e. to search for a solution predominantly using the 
He i diagnostic, we ran test fits with increasing relative weight 
of He i /1447 1 . We found that a solution with an improved fit 
could be obtained for an increase of the relative line weight 
with a factor of five. This solution has a T e g lower by 3.7 kK 
and all other fit parameters approximately equal to the solution 
with the low He i /1447 1 weight. The overal fit quality of the 
other lines has somewhat decreased. As these lines are rela- 
tively strong, this is difficult to discern by eye. Also note that 
within the lower error estimate for r e ff of 3.9 kK the two solu- 
tions agree. Similar results were also obtained for the 02 stars 
BI 237 and BI 253, where a reduction of T e g by 3.8 kK and 
5.4 kK was found, respectively, for an increase in the relative 
weight of He i AAAll by a factor of six and two, respectively. 
Consequently, based on our analysis we can only give a range 
for the r e ff scale of the 02 stars of 49 - 54 kK, the bounds 
being set by the temperatures based on the He i A441 1 and hy- 
drogen/He ii diagnostics. 



The upper end of our hydrogen/helium based 02 effec- 
tive temperature range is in good agreement compared to 
the average effective temperature of 54 kK as determined by 
IWalborn et al] d2004l) . As we already mentioned the giant Nl 1- 
03 1 with a T e g lower by approximately 8 kK compared to the 
other 02 stars forms an exception to this. Its effective tempera- 
ture compares better to the temperature of the 03 star Nl 1-060. 
A close inspection of the spectra of the two stars reveals why 
this is so. Both stars have exactly the same equivalent width 
ratio of the Hei /14471 and Hen /14541 lines, which results in 
similar values for r e ff. To also test whether the relatively strong 
He i /14471 line could be dominating the fit, forcing a relatively 
low r e jf, we refitted the spectrum of Nl 1-031 ignoring the neu- 
tral helium lines. This again resulted in an effective temperature 
of45kK. 



Interestingly, the nitrogen line analysis of IWalborn et al 



for Nl 1-031 did result in a higher effective temperature of 
55 kK. We are not sure whether this discrepancy is the re- 
sult of a systematic offset between the Nv /14603-20/Niv and 
Hei /14471/Hen temperature scales. A more recent compa ri- 
son by one of us (P.A.C.) of the spectral fit of Walb onTet al.l to 
new VLT data, however, showed relatively large discrepancies 
in the helium line fits compared to the nitrogen line fits, indi- 
cating that this could be the case. The discrepancy could also 
be due to diff erences in fitting assumptions and approaches. 



Wal born et al.l adopted a surface gravity of \ogg = 4.0 and 
estimated a mass loss rate of M = 1.0 x lO~ 6 M0yr~' from 
the He n /14686 line. Our fit indicates that the former param- 
eter should be lower by 0.15 dex, consequently lowering 7^. 
With respect to the mass loss rate we find a value higher by 
a factor four. An increase in M of this magnitude can have a 
significant effect on th e strength of different nitrogen lines (e.g. 
ICrowther et al. 2002 ) and, therefore, on the derived effective 
temperature. 



Massev et al.l ([2005 ) analysed a total of 11 02 stars. They 



could determine the effective temperature for three dwarfs and 
one giant, with 47.0 kK < T e g < 54.5 kK. For the supergiants 
only lower limits of T e g > 42 kK and higher could be derived. 
These results are in agreement with our findings. However, 
these authors note that the correlation with T e ff for the 02-3.5 
spectral types is not tight. In particular no good agreement was 
found between the ratios of the N m and N iv emission lines 
and the Hei and Hen lines. Consequently, a more thorough 
investigation of the 02 stars based on both the nitrogen and 
helium spectrum is necessary to resolve this degenerate class 
adequately. 

4.3. Gravities 

The distribution of our programme stars in the log T e g - log g c 
plane is presented in Fig. [2] To calculate the surface gravity 
corrected for centrifuga l accel erati on (logg^) the method dis- 
cussed bv lHerrero et al. (1992) and iRepolust et al.l (2004) was 
adopted. Different luminosity classes are denoted using circles, 
triangles and squares for, respectively, class V, III and I-II ob- 
jects. In this figure we see that the dwarfs, with exception of 
two objects, form a group clearly separated from the latter two 
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Fig. 2. The log T e g - log g c plane for the analysed LMC objects. 
Different luminosity classes are denoted using circles, triangles 
and squares for dwarfs, giants, and bright giants and super- 
giants, respectively. With exception of two objects the dwarfs 
are clearly separated from the luminosity class I-II objects. The 
giants seem to overlap with both the luminosity class V and I-II 
objects. Shown as a dotted line is the average log r e (j - log g c 
relati on of th e SMC O-type stars of luminosity class I-II-III 
fromlPaper 11 



groups. In contrast to our findings in Paper II we do not find 
a clear separation between the dwarfs and giants. Instead the 
latter group of objects shows an overlap with both luminosity 
class V and I-II objects. 

Shown in Fig. [2] as a dotted line is the average log r e jf - 
log g c rela tion of th e SMC O-type stars of luminosity class I- 
II-III from lPaper 11 The majority of the evolved LMC objects 
seem to agree with this trend, that illustrates the evolutionary 
correlation between effective temperature and surface gravity. 
However, note that some evolved objects are at considerable 
distance from the average relation. Thus, a calibration of the 
two parameters fo r a given luminosity class should be taken 
with care (see also lRepolust et al ]|2004 . 

The comparison of spectroscopically determined masses 
(M s ) with masses obtained from evolutionary calculations is 
presented in Fig. [3] Using the same symbols as in Fig.[2]dwarfs, 
giant and supergiants are distinguished. Evolutionary masses 
(M ev ) were derived from the evol utionary tracks calcu lated for 
a metallicity of Z = 0.4 Zq from lSchaerer et aD dl993h . The er- 
rors on these masses correspond to the maximum mass interval 
in the error box spanned by the uncertainties in luminosity and 
effective temperature. 



The tracks from ISchaerer et al.l do not include the effects 
of rotation. Consequently, this additional source of error is not 
included. Calculations including v m show that in some cases 
one can no longer assign an unambiguous M(L, T e ff). This is 
a result of rotationally enhanced mixing or the unknown incli- 
nation angle, if the star has a non-spherically symmetric dis- 
tribution of R+ and T e ff, causing complicate d tracks includin 



loops during the secular redward evolution (IMeynet & Maede 
120051) . As sessing th e impact of u rot on the M ev determination we 
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Comparison of spectroscopic masses with ma s ses de - 

~ dl993l) . 



rived from the evolutionary tracks of ISchaerer et al. 



The one-to-one correlation between the two mass scales is 
given by the dashed line. Symbols have the same meaning as 
in Fig. [2] The objects with the highest evolutionary masses ex- 
hibiting a mass discrepancy correspond to the four 02 stars in 
our sample. 



non rotating tracks to those obtained from tracks calculated for 
y ro t = 300 km s _1 that the error in M ev will not increase by more 
than approximately ten percent. 

In Fig. [3] we see that the majority of the objects are lo- 
cated left of the one-to-one correlation, given by the dashed 
line. The error bars of twelve objects do not even touch 
this correlation. Consequently, we find a significant mass dis- 
crepancy. Simi lar mass discrepancy problems have been dis - 
cuss ed bv e.g. iHerrero et al.l dl992h , (de Koter & Vinkl d2003l) . 
and iRepolust et al.l d2004l) . Most of these classical problems 
wer e attributed to limitations in the stellar atmosphere mod- 
els dHerrerolll993l) and to potential biases in the fitting process 
fsee lPaperB . Here we cannot explain the found discrepancy in 
such a manner. We will provide a more thorough investigation 
and discussion in Sect. [5] 



showed in Paper II from a comparison of masses derived from 



4.4. Helium abundances 

For the SMC sample that we analysed in Paper III we found a 
correlation between the helium surface abundance and surface 
gravity. This could partly be explained by evolutionary effects. 
As the surface gravity decreases when a star evolves away from 
the ZAMS, objects with lower gravities would correspond to 
more evolved objects and are more likely to have atmospheres 
enriched with helium. 

To investigate whether the scenario discussed above also 
applies to the current LMC sample, we plot the helium abun- 
dance as determined with the automated method as a function 
of log g c in Fig. [4] Also shown as a dashed line is a measure 
of the initial helium abundance. This value of Y He - 0.09 was 
calculated by averaging the surface helium abundances of the 
dwarf type objects with a helium abundance smaller than the 
total sample average. Compared to this measure for the ini- 
tial helium abundance a correlation between the surface gravity 
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Fig. 4. Helium abundances as a function of surface gravity. 
Symbols have the same meaning as in Fig. [2] The dashed line 
at iHe = 0.09 is a measure for the initial helium abundance 
and corresponds to the mean of the helium abundances of the 
dwarfs with F H e smaller than the sample average. Starting at 
the highest gravities an increase of the helium fraction is seen 
down to \ogg c =s 3.6. 



and helium enrichment can be observed. Starting at the high- 
est gravities, we find an increase in the average helium abun- 
dance towards lower \ogg c . Note that the two objects exhibit- 
ing the largest helium fractions are supergiants. The increase 
can agai n be partially explained as an evolutionary effect. In 
Paperll a similar correlation between average helium abun- 
dance and gravity was found to exist down to the lowest gravi- 
ties investigated. 

Interestingly, in Fig. H we see that for our LMC sample no 
helium enrichment is found at log g < 3.6. Not even the super- 
giants show evidence of enrichment below this gravity. Why is 
this so? In Fig. [5] where we present the HR-diagram for our 
sample, the answer to this question is given. Shown as a grey 
area in this figure is the region in which the rotating evolu- 
tionary models of iMeynet & Maederi (12005b predict a helium 
surface enrichment of at least ten percent. The majority of the 
evolved objects are located outside this region. Consequently, 
based on their specific evolutionary phase these objects are not 
expected to show any enrichment. Also important is the fact 
that we have selected the hottest objects from the Nil cluster. 
As a result of this our sample is biased and does not contain 
the low gravity objects with a high luminosity, i.e. those more 
likely to be enriched, as these have evolved into cool B-type 
stars. 

Apart from the supergiants also four dwarfs are found to be 
enriched. Figure |3J in which we have highlighted stars with a 
helium abundance of at least 0.12 using open symbols, shows 
that one these dwarfs, Sk -66 18, is located relatively close to 
the region in which enrichment is predicted. Consequently, ro- 
tationally enhanced mixing is a possible explanation. The three 
remaining dwarfs, Nl 1-60, Nl 1-065 and Sk -70 69, in contrast 
lie relatively close to the ZAMS. Therefore, "normal" mixing 
cannot explain their enrichment. Instead, a possible explana- 
tion is given by chemically homogeneous evolution. As this 




Fig. 5. Hertzsprung-Russell diagram for the LMC sample. 
Symbols have the same meaning as in Fig. [2] Over plotted as 
grey lines are the evolutionary tracks of ISchaerer et al. ( 1993 ) 
for Z = 0.4 Zq, with a black line representing the ZAMS. Open 
symbols indicate objects with F H e ^ 0.12. The grey area corre- 
sp ond to the region in whic h the rotating evolutionary models 
of IMeynet & Maederi d2005l) predict a relative helium enhance- 
ment of at least 0.01. 



may also be linked to the mass discrepancy, we will return to it 
in Sect.|5] 



4.5. Microturbulence 

The microturbulent velocities determined with the automated 
fitting method are also given in Tab. [2] Though the error es- 
timates are relatively large (see Tab. [3]), we find that for the 
current data set the y tur b determinations were sufficiently accu- 
rate to reveal a weak correlation between this parameter and 
the surface gravity. This is shown in Fig. [6j where it can be 
seen that for log g < 3.6 the average microturbulence recovered 
from the line profiles increases systematically. The situation for 
log g > 3.6 is less clear, as the error bars are on average larger 
and the values for are more or less randomly distributed 
between and 20 km s" 1 . The reason for this is that for larger 
values of the surface gravity the line profiles become intrinsi- 
cally broader due to the increased Stark broadening, making it 
more difficult to accurately recover i> tU rb from the line profiles 
only. 

Based on samples predominantly consisting of unevolved 
early B-type Galactic stars other authors have also found 
a relation b e tween microturbulence and su rface gravity, e.g. 
Kili an et al.1 dl99ll) . iGies & Lambert! d 19921) and iDaflon et al 



d2004 . More recently Hun ter et alj 12006) analysed a sample 
of early B-type stars in the Magellanic clouds and also found 
a trend of increasing y tur b for decreasing log g. To derive the 
values for u tur b all these authors relied on curve-of-growth tech- 
niques, which were applied to metal lines such as those of Si m 
and O in calculated using plane parallel models. Consequently, 
our line profile based analysis is an independent confirmation 
of the existence {or requirement by lack of a physical explana- 
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Fig. 6. Microturbulent velocities determined using line profile 
fits as a function of surface gravity. Symbols have the same 
meaning as in Fig. [2] For loggr < 3.6 a trend of increasing 
Uturb with decreasing log g is visible. For larger gravities the 
uncertainties in the v tur \, determinations are too large to discern 
any possible relation. 
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Fig. 7. Microturbulent velocity as a function of the location of 
the line forming region of Hei AAAll, which is defined as the 
location where the radial optical depth in the line core reaches 
a value of r = 2/3. A weak trend is visible that suggests that 
for increasing extension of the atmosphere larger values of u tur b 
are necessary to fit the line profiles. 



Hon or failures in the line broadening mechanisms) of micro- 
turbulence in the atmospheres of these type of stars. 

The physical mechanism explaining the observed microtur- 
bulence or its relation to t he s urface gravity is stil l poorl y un- 
derstood. |Kudritzkj ([1992]) and lLamers & Achmadld 19941) have 
argued that the observed microturbulence might be the result of 
a stellar outflow, implying the existence of a velocity gradient 
in the photospheric layers , which can mimic micro turbulence- 
like desaturation effects. ISmith & Howarfhl d 1998b . however, 
showed by applying a simple core-halo model to the Galactic 
09.7 supergiant HD 152003 that this effect would not be suffi- 
cient to explain the observed u tUT \,. Indeed, our analysis employ- 
ing a unified photosphere and wind model confirms that the mi- 
croturbulence cannot be explained as an artifact of a transonic 
velocity field. 

A possible explanation for the gravity dependence could 
be related to instabilities in the wind. These instabilities are 
reflected in the large turbulent velocities (~100-200 kms -1 ) 
necessary t o fit the wind lines in the inte r mediate and outer 
wind (e.g. Groene wegen & Lamers 1989: lHaseretai]ll998l 



Evans et al. 2004b) and could be related to shocks due to the in- 



trinsic line-driven ins tability dLucvll 1983c lOwocki et al.M 1988 



Owocki & Puisll999t) . For a low surface gravity the wind starts 



at larger Rosseland optical depth compared to the high grav- 
ity case. Therefore, the line forming region of low gravity ob- 
jects contains a relatively large contribution originating from 
(the base of) the wind. Consequently, this region could be af- 
fected by the onset of the line-driven instability introducing 
wind-turbulence into the line profiles. 

Tentative support for the above described scenario may be 
implied by Fig. [7] where we show o tur b as function of the line 
forming region of He i /14471 in units of the stellar radius. The 
location of the line forming region is defined as the position 
at which the radial optical depth in the line core reaches a 
value of 2/3. This figure shows that when the radial distance 



to this position increases also the average microturbulent ve- 
locity increases. Though the trend is weak, it does seem to 
indicate that when the atmosphere becomes more extended, 
higher values of i> tur b are necessary to reproduce the line pro- 
files. It may appear that this implies that the line forming re- 
gion enters in to the regime where wind turbulence develops. 
However, for this statement no compelling evidence is avail- 
able, as we do not find any correlation between u tul -b and the 
distance between the line forming region and (for instance) the 
sonic point. Moreover, for most objects at R{T\i ne = 2/3) > 1 .04 
values of u tur b close to the maximum allowed value are found. 
Consequently, these values should be interpreted as lower lim- 
its. Therefore, we can only conclude that our analysis points 
towards a gradient in the turbulent velocity, possibly connected 
to a link between microturbulence and wind instabilities, and 
suggest further investigation in this direction. 

4.6. Wind parameters 

The wind parameters and their uncertainties determined using 
the automated method are listed in Tabs. [2] and [3] Compared 
to our SMC analysis we find that we were able to accurately 
determine these parameters for a significantly larger number 
of objects. This is mainly the result of an on average higher 
signal-to-noise ratio of the spectra as well as of denser winds 
for the LMC objects compared to their SMC counterparts. In 
total we determined 22 mass loss rates and 6 upper limits. The 
upper limits are defined (and can be identified in Tab. [3]) by an 
error bar - log M > 1 .0 dex. 

To provide a meaningful comparison of the mass loss rates 
we place the LMC objects in the modified wind momentum 
luminosity diagram. This diagram shows as function of stellar 
luminosity the distribution of the so-called modified wind mo- 

1 /2 

mentum, which is defined as D mom = Mv^RJ . Not only does 
this allow for an assessment of the behaviour of M within our 
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Fig. 8. Modified wind momentum (D mom ) in units of 
gcms^SQ 1 ' 2 vs. luminosity. Symbols have the same mean- 
ing as in Fig. [2] Upper limits are shown as inverted triangles. 
Dashed lines correspond to the p redicted w i nd-momentu m lu- 
minosity relations (WLR) from I Vink etaD d200d l200lh . The 
upper, middle and lower relation, respectively, correspond to 
predictions for Galactic, LMC and SMC metallicity. The ob- 
served modified wind momenta show a strong correlation with 
luminosity with an average relation that lies in between the pre- 
dicted Galactic and SMC WLR. This is quantified by the em- 
pirical WLR that was constructed for the LMC objects and that 
is shown by the dotted line. The open square corresponds to 
the wind momentum corrected for clumping of the supergiant 
Sk -67 166. Shown as a dashed-dotted is the empirical WLR 
obtained using this corrected D mom . 

sample, it also provides a convenient method to compare the 
observed wind strengths to the predictions of line driven wind 
theory. According to this theory D mom is predicted to behave as 



log D,r 
where 



xlog (L*/L©) + log D 



(1) 



Puis et al. 



is t he stellar luminosity ( Kudri tzki et al.l 1995 



1996). In this equation x is the inverse of the slope of 



the line -strength distribu tion function corrected for ionisation 
effects dPuls et alJ hoOO). The vertical offset D a is a measure 
for the effective number of lines contributing to the accelera- 
tion of the outflow. 

In Fig. [8] the distribution of the modified wind-momenta 
for our programme stars are presented. Indicated using circles, 
triangles and squares are objects of, respectively, luminosity 
class V, III and I-H. Upper limits are shown as grey inverted 
triangles. A clear correlation between L* and D mom can be ob- 
served in this figure. Over an order of magnitude in the av- 
erage modified wind-momentum decreases by approximately 
1.5 dex. A comparison of the behaviour of D mom with pre- 
dictions is facilitated b y the theoretical WLRs calculated by 
Vink et al. I d2000l 2001) that are shown as a set of dashed lines. 
The upper, middle and lower of these predicted power laws 
were calculated for, respectively, Solar, LMC and SMC metal- 
licity. Compared to these predictions we find that the LMC 
wind-momenta are approximately confined between the theo- 
retical WLR for Galactic and SMC metallicity. In other words, 



compared to the Galactic and SMC case, stars in the LMC have 
intermediate wind strengths. 

To further quantify the behaviour of LMC winds relative to 
that of Galactic and SMC outflows we have fitted a power law 
to the observed modified wind-momentum distribution, consis- 
tently accounting for both the symmetric errors in L* and the 
asymmetric errors in D mom , This yielded the following empiri- 
cal WLR 



logD mom = (1.81 + 0.18) log (L*/L G ) + (18.67 ± 1.01) . (2) 



The theoretical LMC relation from Vin k et al.l is given by 
x = 1.83 and D mom = 18.43. The error bars of the theoretical 
and empirical relations are in agreement. More importantly, in 
Fig.[8]the empirical relation, shown as a dotted line, is found to 
lie between the predicted Galactic and SMC relations. These 
predictions have been foun d to be in good agreement with 
the observed Galactic WLR ( Repolust et al . 2004, iPaper J) and 
observed SMC WLR dPaper III) . Consequently, our empirical 
LMC WLR is quantitative evidence for the fact that massive 
stars in this system have mass loss rates intermediate between 
those of massive stars in the Galaxy and SMC. 

The differences between the empirical and theoretical LMC 
WLR at the start and end of the observed luminosity range 
are, respectively, 0.17 and 0.16 dex. In Fig. [8] this seems to 
imply a systematic offset between the two relations. However, 
we note that these differences are still smaller than the typi- 
cal uncertainty in D mom of 0.2 dex. More importantly, no cor- 
rection was applied for the possibility that the winds of our 
sample stars are, in contrast to our assumption, not smooth but 
structured. In recent years evidence has been mounting for this 
so-called clumping in the winds of O- and early B-type stars. 



( Crowther et al. 2002 Iffillier et alJ 12003 


: BouretetalJ 2003: 


Massa et alJ 120031: iMartins et all [2004, 


2005bt Bouretetal. 


20051: iFullerton et al]l2006h seems to suggest the existence of 



clumping factors in the range of 10-100, implying correspond- 
ing reduc t ions o f M by factors 3 up to 10. More recently, 
IPuls et alJ d2006l) also found clumping factors of the order 5 



to 10 (normalized to the unknown clumping properties in the 
outermost, radio-emitting wind) from the analysis of Ha, in- 
frared, millimetre and radio fluxes. In the present study we try 
to account for possible wind clumping effects b y correcting the 



mass loss rates of stars with Ha in emission. iMarkova et al 



( 2004) and Re polust et alJ d2004h argue that the mass loss rates 
of these stars could be overestimated as a result of the fact that 
Ha emission lines are formed over a relatively large volume 
where clumping might have set in. In contrast, for stars with 
Ha in absorption the line is formed relatively close to the stel- 
lar surface, where clumping effects are negligible. Based on the 
comparison of dwarfs and supergiants in their Galactic sample 
iRepolust et al. derived a numerical correction factor of 0.44 for 
the mass loss rates of supergiants with Ha in emission. 

We have applied the clumping correction to the super giant 
Sk -67 166, which has a Ha emission profile. In Fig.[8]its new 
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wind momentum is indicated using an open symbol. Using this 
value the following empirical WLR is obtainecQ 

logD mom = (1.43 ± 0.17) log (L*/Z©) + (20.77 ± 0.97) . (3) 

In Fig. [8] we see th at for log L + / > 5.3 the new WLR com- 
pares better to the Vink et all relation. For lower luminosities 
the situation is less clear. Due to the large uncertainties, this 
range has a relatively low weight in the fit. Consequently, a 
discrepancy for low L+ is less significant than the good agree- 
ment obtained for the higher luminosities. For this reason it 
is difficult to investigate the existence of a "wea k wind prob- 
lem" f or stars at log L+/Lq < 5.3, first reported by B ouret et all 
d2003l). This study as well as later studies dffillier et al.ll2003t 
Evans et al]l2004at iMartins et al.ll2004 l2005bl) report a steep- 
ening in the WLR relation relative to predictions starting at 
about the above mentioned L+, leading to an over prediction of 
the wind strength by up to a factor 100 at logL*/L© ~ 4.5. 
Our LMC results do not appear to confirm this break be- 
tween observations and seem to follow the predictions down 
to log L+ILq =s 5.0. In a forthcoming paper we will present a 
comprehensive overview of the observed WLR relations in our 
Galaxy and the Magellanic Clouds as well as a thorough dis- 
cussion of the successes and failures of the theory of radiation 
driven winds in predicting the WLR, including possible causes 
for the weak wind problem (Mokiem et al. in preparation). 



5. The mass discrepancy 

In Fig. [9] we investigate the mass discrepancy of our LMC tar- 
gets as a function of the helium surface abundance. On the ver- 
tical axis a measure for this discrepancy is shown, which is 
scaled to the mean of the evolutionary and spectroscopic mass. 
This ensures that positive and negative discrepancies follow 
the same linear scale. For non-enriched stars, i.e. Yu e < 0.10, 
approximately three times as many objects lie above the one- 
to-one correlation than below it. Therefore, in contrast to our 



finding for the equivalent SMC case ( {Paper I A) a significant 
mass discrepancy is found for our sample of non-enriched LMC 
stars. The reason for this is unclear. As we stated before our 
analysis employs state-of-the-art atmosphere models, and is in 
principle not hampered, as were previous studies, by poten- 
tially unoptimised fits. As stated, in our SMC dataset, analysed 
in an identical manner, no evidence was found for a mass dis- 
cre pancy for Yu e < 0.10 . 

Masse v et al. I (120051) also study the M s vs. M ev problem in 
a set of LMC stars. For objects hotter than 45 kK they find 
a mass discrepancy that is even stronger than what we find. 
This behaviour could be the result of an underestimate of the 
photospheric line pressure in this high temperature regime. In 
Fig.|9]we have highlighted the objects with T^ff > 45 kK using 
open symbols. Although their mass discrepancy is considerable 



1 The small errors in the relevant parameters of Sk -67 166, which 
dominates the WLR at very high luminosity, cause the significant dif- 
ference between Eqs. 2 and 3. Note that our least square fitting does 
not account for the correlation of both quantities (due to R+). In view 
of the well known dis tance these effects are likely small compa red to 
the Galactic case (see lMarkova et al.l2 004; Rep olust et"ai]|2004l) . 
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Fig. 9. Mass discrepancy as a function of helium abundance 
for the LMC sample. Evolutionary masses used to calculate 
the discrepancy were derived from the non-rotating tracks of 
ISchaerer et al.l dl993l) . Symbols have the same meaning as in 
Fig. [2] The open symbols correspond to stars with r e ff > 45 kK. 



they do not stand out as a separate group. To further illustrate 
this, note that the hottest supergiant in our sample Sk -67 166 
(r eff = 40 kK) at Y He = 0.28 has a M s of 75 M© which is 
in very good agreement with its M ev of 70 Mq. Consequently, 
though we can not explain the mass discrepancy we do not an- 
ticipate that it is connected to a flawed treatment of the photo- 
spheric radiation pressure that manifests itself at metallicities 
as high as that of the LMC environment (but not yet at values 
typical for the SMC). 

If we accept the analysis of stellar and photospheric pa- 
rameters, the presence of a mass discrepancy may point to an 
oversimplified picture of the evolution of massive stars used 
to determine M ev and/or, possibly, to a breakdown of the as- 
sumption of a spherically symmetric atmosphere. One obvious 
simplification may be that we have used tracks for non-rotating 
stars. One of the effects introduce d by rotation i s a wide bi- 
furcation in the evolutionary tracks (Maed erll987 ). Stars rotat- 
ing faster than roughly half the surface break-up velocity will 
essenti ally follow tra cks representative of homogeneous evo- 
lution. iLangerl d 19921) showed that as a result of this the M/L- 
ratio will be a monotonically decreasing function for increasing 
helium enrichment. Consequently, stars evolving along homo- 
geneous tracks are expected to be increasingly under-massive 
for increasing helium abundance. One would therefore derive a 
positive mass discrepancy if the evolution of a rapidly rotating 
star was incorrectly described using non-rotating or modestly 
rotating model tracks. 

So, can rapid rotation be used to explain the observed mass 
discrepancy? Let us first focus on the supergiants. For these 
stars the mass discrepancy problem appears absent, though we 
note two exceptions at Fne = 0.07 and Yu e = 0.09. This ab- 
sence of signs of a distinct mass discrepancy was also foun d 
in our study of SMC bright giants and supergiants (Pap er lib . 
The conclusion seems to be that the supergiants follow non- 
rotating or modestly rotating evolutionary tracks. This implies 
that, apparently, all supergiant targets that we have selected in 
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both the SMC (IPaper II) and LMC happen to have started out 
their evolution with low or modest initial rotational velocities. 



As the initial rotational velocity distribution derived in Paper II 



for the NGC 346 cluster implies that only some 5-15 percent 
of stars are expected to evolve along homogeneous tracks this 
should not be alarming. One supergiant, however, might not fit 
this picture: Sk -67 166. This source does not feature a mass 
discrepancy, but does show an increased helium abundance. 
Given its position in the HRD this is quite difficult to explain in 
term of evolution without rotation. A scenario accounting for a 
more delicate interplay between mass loss and rotation may b e 
required to explain its properties (see lHerrero & Lennonl20 04). 

Now let us turn to the dwarfs and giants. For our SMC 
sample we found a correlation between the helium surface en- 
richment and the mass discrepancy for class V and III stars, 
with M s being systematically smaller than M ev for Fee 0.11. 
Interestingly, this dependence was only found for dwarfs and 
giants (see above). It was suggested that this behaviour was 
the result of rotationally enhanced mixing, enriching the at- 
mospheres with primary helium. The dwarfs and giants at 
Yne > 0.11 displayed in Fig. [9] also systematically suffer from 
a positive mass discrepancy. Consequently, this scenario also 
seems to be a good explanation for the situation in the LMC 
sample. However, we realise that these seven objects are a rel- 
atively small sample and that the typical uncertainties in Fn e 
are 0.03. This statement should therefore be seen as a working 
hypothesis. 

For helium enriched dwarfs and giants the above scenario 
seems a logical one, as helium enrichment early on in the evo- 
lution must imply efficient mixing. For the dwarfs and giants 
that do not show excess helium in their surface layers such 
a clue or indication is not present (unless they evolve left of 
the main sequence as do chemically homogeneous stars, see 
Sect. l5.11 >. In principle these non Fee-enriched stars could ro- 
tate sufficiently rapid to cause a mass discrepancy, but not to 
the extent of chemically homogeneous evolution. However, if 
so, a very significant fraction of the stars should have started 
their evolution at super-critical rotatio n. At lea st for the stars 
studied in the SMC cluster NGC 346 dPaper lib this is not the 
case. 



5.1. Chemically homogeneously evolving 02 stars 

Three out of the four 02 stars are found to lie to the left of the 
ZAMS in Fig. [5] and exhibit relatively large mass discrepan- 
cies (see open symbols in Fig.|9]l. These are Nl 1-026, BI 237 
and BI 253. For a discussion of the effective temperature of the 
fourth 02 object, Nl 1-031 (r eff ~ 45 kK), we refer to Sect.gj] 
Even though the hottest three (r e (j ~ 53 - 54 kK) are not sig- 
nificantly enriched in helium it is tempting to speculate on a 
possible (near) homogeneous evolution of these objects. This 
would not only provide an explanation for their mass discrep- 
ancy, it would also explain their peculiar location in the HRD. 
The possibility that these are true ZAMS stars is very exciting. 
However, due to the short evolutionary time scales in this part 
of the HR-diagram, it w ould also be unl i kely . Based on the non- 
rotating model tracks of lSchaerer et al.1 d 1993b we estimate that 



a 60 Mq star will within 2 Myr evolve away from the ZAMS 
to a location in the HRD at T e ff » 45 kK. This is well beyond 
the error bars on the parameters of these objects. Given the fact 
that Nl 1-026 and BI 253 are thought to be associated with the 
LH 10 and 30 Doradus clusters, respectively, w hich have ages 
of ap proximately 3 (see Sect. 16.21 ) and 2 Myr Jde Koter et al 



Il998l) . a normal evolutionary scenario appears unlikely. Based 
on its large r adial velocity the fie ld star BI 237 probably is a 
runaway star (iMassev et al. 2005 ). therefore it is also likely to 
be relatively old. 

The problems with a reconciliation of the hottest three 02 
stars with fully homogeneous evolution and an age of at least 
2 Myr, are i) that their surface helium abundances are not sig- 
nificantly enriched, and ii) that their rotational velocities are 
not extreme. The latter issue need not be "a smoking gun" con- 
sidering the possibility that we may see them relatively pole- 
on and the fact that their u r sin i values, ranging between 110- 
190 km s -1 , are above the sample average. Concerning the first 
point, it appears that we must concede to the possibility that 
stars may evolve along tracks similar to those for homogeneous 
evolution while in fact they are not fully homogeneous - i.e. 
the near surface layers are not yet strongly affected by the mix- 
ing that must occur deeper in. Having said this, we do point to 
the fact that the error bars on Fn e do allow for relatively large 
ages even within the hypothesis of fully homogeneous evolu- 
tion. Within the error bars the maximum Fn e is ~ 0.13 for all 
three stars. It takes fully homogeneously evolving stars of 60 
and 40 Mq about ~ 1.5 and ~2.0 Myr, respectively, to build up 
this amount of helium enrichment (S. -C. Yoon, private com- 
munication). This is close to the derived cluster ages of LH 10 
and 30 Doradus. For reference, if the star would evolve along 
near-homogeneous tracks they would be older. 

We tentatively conclude that the HRD position, helium 
abundance, and rotational velocity of the three hottest 02 stars 
in our sample could be consistent with (near-)homogeneous 
evolution. Evidence for efficient rotation-induced mixing dur- 
ing t he main sequence p hase of O stars has also been presented 
bv lLamers et alj (1200 lb . based on the chemical abundance pat- 



tern of ejected circumstellar nebulae. 
5.2. Large sample trends 

To firmly establish our findings with respect to the mass dis- 
crepancy problems, we compare in Fig. [10] spectroscopic and 
evolutionary masses for all stars that have been analysed using 
our automated fitting method so far. Pres ented in this figure are 
the combined Galactic (Paper I), SMC (IPaper lib and the cur- 
rent LMC samples, corresponding to a total of 71 O- and early 
B-type stars. 

The figure confirms our two main findings. First, the good 
agreement between M s and M ev for the supergiants is clearly 
visible. This is an encouraging finding and we believe that this 
means that the improvements in the stellar atmosphere mod- 
els, evolutionary calculations and spectral analysis techniques 
have fina lly resolved the lon g standing mass discrepancy as 
found by iHerrero et al.1 (11992b . The correspondence between 
the spectroscopic and evolutionary mass scale, in particular for 
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Fig. 11. The FLAMES field for the star forming region Nil with the associations LH9 (south of centre) and LH1 (north of 
centre). The FLAMES targets are identified using circles and their star identification is given. From lEvans et al. (2006). [Included 
as separate JPG file in astro-ph submission.] 
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Fig. 10. Mass discrepancy as a fu nction o f helium ab undance 
for the combined samples from iPaper II (Galactic), iPaper II 
(SMC) and the current LMC sample. Symbols have the same 
meaning as in Fig. [2] Note the good correspondence between 
M s and M ev for the bright giants and supergiants (square sym- 
bols). 



increasing helium enrichment is striking. The reason for this is 
probably related to the fact that the class I-II objects are found 
in the region of the HR-diagram in which stars are not expected 
to undergo extreme evolutionary phases, such as homogeneous 
evolution. Consequently, the enriched bright giants and super- 
giants have evolved along relatively simple evolutionary tracks, 
which (in this respect) appear well understood. 

Our second finding, i.e. that of a correlation of mass dis- 
crepancy with the helium abundance, is also corroborated by 
the large sample in Fig. [TO] For F He > 0.10 all dwarfs and gi- 
ants, except for one dwarf, are found above the M s - M ev line. 
In all fairness, given the typical uncertainty of 0.03 in Yu e , from 
a statistical point of view the region 0.09 < Yu e < 0.12 should 
be regarded with care. For larger helium abundance, however, 
the correlation can be regarded to be statistically significant. A 
total of 1 1 objects show a positive mass discrepancy, with only 
a single counter example. Moreover, the magnitude of the dis- 
crepancy seems to be related to the amount of enrichment. This 
strongly points to efficient mixing in the main sequence phase, 
leading to (near-)chemically homogeneous evolution. 



6. The evolutionary status of N1 1 

In this section we explore the evolutionary status of the Nil 
field. We will first briefly outline the current understanding of 
Nil with respect to the OB associations LH9 and LH10 and 
in particular the sequential evolutionary link between the two. 
Based on the analysis of the 22 stars in our sample associated 
with these clusters we will then estimate their ages and discuss 



whether they are compatible with a sequential star formation 
scenario. 



6.1. LH9andLH10 inN11 

Nil is an intricate giant Hn region containing several massive 
star forming complexe s. The largest of these a re the OB associ- 
ations LH9 and LH10 jLucke & Hodgell970h . An image of the 
region, inc luding the adopted star identifications, is presented 
in Fig. [TT1 dEvans et al.ll2006l) . Several studies have suggested 
that the star forming activit y in these a ssociations are linked. 
In particular the study by iParker et al. I dl992l) has provided a 
key to understanding of the structure and formation of LH9 
and LH10. Their analysis of the stellar content of Nl 1 revealed 
the presence of several 03-05 stars and possible ZAMS stars 
in LH10. In contrast the earliest spectral type associated with 
LH9 was found to be 06. They also found the slope of the ini- 
tial mass function of LH10 to be significantly flatter than that 
of LH9, indicating that the former contains a higher ratio of 
high mass to low mass stars. Combined with the fact that the 
reddening of LH10 is larger than that of LH9, it was concluded 
that LH10 is the younger of the two clusters. Based on these 
findings the authors propose an evolutionary link between the 
two, where star formation in LH10 could possibly be triggered 
by the stellar winds and supernovae of massive stars in LH9. 



Further evidenc e for a sequentia l star formation sce- 
nario was given by IWalborn & Parked (119921) . who found a 
dual structural morphology of Nil analogous to that of the 
30 Doradus H n region. In the latter substantial evidence sug- 
gests the presence of current star formation in regions sur- 



rounding th e centra l star cl uster (e.g . Walborn & Bladesl ll987 
Hyland et all 19921 also see|Walborn & Bladesll997l) . This sec- 



ondary burst seems to be set off approximately 2 Myr after the 
initial star formation took place, possibly i nitiated by the en- 
ergetic activity of the evolving cluster core. IWalborn & Parker 
argue that this is very similar to Nl 1 where the stellar content of 
LH 9 and LH10 also sugg ests an age difference of ~2 Myr (also 



see IWalborn et al.l ll999). The process in Nil, however, would 
be advanced by ~2 Myr, classifying it as an evolved 30 Doradus 
analogue, though less massive. 

More recently, several bright IR sources showing charac- 
teristics of young stellar obj ects were discovere d in the N11B 
nebula surrounding LHlOby Barba et al.l (120031) . These objects 
are probably intermediate mass (pre-)main-sequence Herbig 
Ae/Be stars belonging to the same generation as do the LH10 
objects (the pre-main-sequence evolutionary timescales of in- 
terme diate mass s tars being longer than that of massive OB 
stars). iBarba et aL also found that the massive stars in LH10 
have blown away their ambient molecular material and are cur- 
rently disrupting the surface of the parental molecular cloud 
material surrounding LH10 (i.e. the material in Nl IB). 
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Fig. 12. Comparison of the programme stars located in N 1 1 with isochrones derived from the evolutionary tracks of lSchaerer et al 
(1993). The symbols used to denote dwarf, giants and bright- and supergiants are, respectively, circles, triangles and squares. 
Filled and open symbols are used to distinguish stars associated with LH9 and LH10, respectively. Black symbols refer to objects 
located within a radius of two arc minutes of the LH9 and LH10 core, whereas the grey symbols correspond to stars found outside 
this radius. Labels correspond to the Nil identification numbers given in Tab. Q] Is ochrones (dashed lin es) are shown from one 
up to ten million years with 1 Myr intervals. For reference evolutionary tracks from lSchaerer et aL ( 1993 ) are drawn as grey lines, 
where for clarity purpose only the blueward part of the evolution for the most massive tracks is shown. The ZAMS corresponding 
to these tracks is given by the black solid line. 



6.2. The ages of LH9 and LH10 



To estimate the ages of our programme stars we compare their 
location in the HR-diagram to theoretical isochrones in Fig. [12] 
The luminosity classes V, III and I-II are represented using, re- 
spectively, circles, triangles and squares. To differentiate be- 
tween stars associated with LH9 and LH10, respectively, filled 
and open symbols are used. Membership is defined on the ba- 
sis of minimum distance to either the LH9 or LH10 cluster 
core. This, therefore, should be taken with some care: the cores 
are about 4 arcmin or 60 parsec apart, which can be traversed 
in 2 Myr if the proper motion of the star is some 30 kms~'. 
Runaway O- and B-type stars have typical velocities of 50- 
100 kms~', therefore, it is entirely possibly that (a few) indi- 
vidual objects are assigned to the wrong association. For posi- 
tion reference see Fig.QT] Isochrones in Fig. [12] are shown as 
dashed lines for one up to ten million years with 1 Myr intervals 



and we re derived from the evolutionary tracks of lSchaerer et al 



( 1993). Note that these tracks do not account for the effects of 
rotation. 



The distribution of the stars in Fig.[T2]is such that they can 
be separated in a group of objects younger than three million 
years and a group of objects older than this age. We also see that 
the oldest objects are predominantly found in LH9, whereas 
LH10 contains the largest fraction of young objects, suggest- 
ing that the clusters can indeed be separated in terms of age. At 
first sight it seems that they can be characterised by an age of 
~7 Myr and ~2 Myr, respectively. However, a significant num- 
ber of objects in both clusters are found to be several million 
years younger and older than these preliminary ages. A possi- 
ble explanation for this large age scatter could be "contamina- 
tion" by field stars. To assess this possibility we have assigned 
grey symbols to all objects outside a radius of two arc minutes 
from the cluster cores. Disregarding these objects reduces the 
age scatter significantly; still a number of stars appear to con- 
tradict with the notion of two coeval populations. 

To investigate the age distributions in more detail the indi- 
vidual age estimates are shown in Fig. Qj] LH9 and LH10 ob- 
jects, which are shown using the identical symbols as in Fig.fTZl 
are placed in, respectively, the left and right part of the diagram 
and are separated using a dashed line. First concentrating on the 
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Fig. 13. Individual ag e estimates for stars in Nl 1 based on the 
non-rotating tracks of lSchaerer et al.l (11993b . Symbols have the 
same meaning as in Fig. [T2] The left and right part of the 
diagram, respectively, contain stars associated with LH9 and 
LH10 and are separated using a dashed line. The upward point- 
ing arrow denotes a lower age limit for Nil -065 determined 
from its surface helium abu ndance using chemically homoge- 
neous evolutionary models dYoon et al. 2006 : Yoon & La nger 
2005). 



LH9 objects we see that they are near-coeval with exception of 
four dwarfs. Of these the two non-core members are located 
at a distance of approximately six arc minutes from the central 
concentration. Therefore, it is probable that they are either spa- 
tially not related to Nil or that possibly their formation was 
triggered more recently by the stellar activity in LH9. The third 
dwarf Nl 1-065 occupies a location very close to the ZAMS and 
in principle would be the youngest member of LH9. However, 
we also find a considerable helium enrichment for this object 
(T He = 0.17). Combined with the fact that Nl 1-065 has a large 
mass discrepancy, this suggests that it might be evolving chem- 
ically homogeneously. Consequently, a more appropriate age 
estimate should be derived from tracks a ppropriate for this kind 



of evolution. Adopting such tracks from lYoon et al. I (120061 also 
see I Yoon & La nger 2005) we derive a lower limit of 6 Myr for 
the age of Nl 1-065 based on the surface helium abundance. 
In Fig. [13] this estimate is indicated using an upward pointing 
arrow, and is in good agreement with the bulk of the LH9 stars. 

The dwarf Nl 1-123 exhibits no chemical peculiarities and 
its location on the sky places it in the central concentration of 
LH9, suggesting that it should have been formed in the burst 
of star formation that formed the cluster. However, given the 
fact that Nl 1-123 is in principle the only star deviating from 
the coeval nature of LH9 we suspect that it is like Nl 1-058 
and Nl 1-066 part of the periphery of the cluster and that its 
relatively close position to the core is due to a projection ef- 
fect. Consequently, we conclude that the central concentration 
of LH9 is coeval with an age of ~7.0 +1.0 Myr. The error bar 
is of the same order of magnitude as the 0.5-1.0 Myr intro- 
duced by uncertainties in the evolu tionary tr acks due to effects 
of (relatively modest) rotation (see lPaper III) . 



The right part of Fig. [13] shows that the stars in the cen- 
tral concentration of LH10 have ages ranging from one up to 
approximately six million years. Despite this large scatter it is 
clear that the majority of the stars are younger than 4.5 Myr 
and that only one object (Nl 1-087) in the cluster core is older 
than this age. An explanation for the large age of the latter ob- 
ject might be that it formed in LH9 and, over time, migrated 
towards LH10. As explained at the start of this section, this 
is a possibility given that the distance from the centre of LH9 
to the current position of Nl 1-087 is only three arc minutes. 
Considering the error bars on the age determinations and the 
addi tional un certainty of 0.5-1.0 Myr introduced by rotation 
(see lPaper It) , we finally estimate an age best describing LH10 
of -3.0 ±1.0 Myr. 

6.3. Sequential star formation? 

If the starbursts in Nil are sequential the formation of LH10 
should be induced by supernova explosions and/or stellar winds 
in LH9. Given the age difference and distance between the 
two clusters we find that both means of triggering are possi- 
ble. In the first case the age difference of approximately four 
million years is compatible with the time of ~3 Myr it takes 
before the most ma ssive stars end their life in a supernova (e.g. 
Schaller et al.ll 1992 ). The time needed for the supernova shock 
to cross the distance of ~3 arc minutes between LH9 and LH10, 
corresponding to approximately 40 par sec at the d istance of 
the LMC, is only ~10 5 years (see e.g. lFalldll981l) : hence it 
is a possible scenario. If we consider trigg ering through stellar 



winds , the time scales are also compatible. iGarcia-Segura et al 
dl996l) . for instance, using hydrodynamical simulations have 
shown that a 60 Mq star can create a wind-driven bubble in the 
interstellar medium of ~50 pc during its main sequence lifetime 
of 3 Myr. Consequently, this scenario seems to be appropriate 
as well and in agreement with the determined time scales. 

In view of the above we conclude that a sequential scenario 
for LH9 and LH10 seems very likely. A combination of super- 
novae and stellar winds from stars in LH9 may have initiated 
star formation in LH10. 



7. Summary and conclusions 

We have analysed a sample of 28 massive OB-type stars lo- 
cated in the LMC. For a homogeneous and consistent treatment 
of the da ta we employed the automated fitting method devel- 
oped bv lMokiem et al. (2006), which com bines the genetic al- 
gorithm based optimisation routine pikaia (Charbonn eaul 1995 ) 
with the fast non -LTE unified stellar atmosphere code fastwind 
( Puis et ap2005 ). The sample is mostly drawn from the targets 
observed with in the context of the VLT-FLAMES survey of 
massive stars (Evans et al. 2005). In total 22 of these stars are 



located in the LH9 and LH10 clusters within the giant Hn re- 
gion Nil. This region is believed to have been the scene of 
sequential star formation, with the stellar activity in LH9 ignit- 
ing secondary starbursts in different associations in Nil. Our 
main findings are summarised below. 
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i) The effective temperature per spectral sub-type of the LMC 
stars is found to be intermediate between that of Galactic 
and SMC O- and early B-type stars, with the LMC objects 
being, respectively, cooler and hotter by typically ~2 kK. 

ii) Based on the helium and hydrogen lines it was possible to 
determine the effective temperatures, though with relatively 
large error bars, of the four 02 stars in our sample. Three 
of these are found to be hotter by more than 3 to 7 kK (see 
Sect. l4.2b compared to the 03 star in our sample, suggesting 
that 02 stars indeed represent a hotter subgroup within the 
O-type class. However, we note that for the one ON2 star 
a relatively low r e jf of 45 kK was obt ained, indicating tha t 
the N iv and N m classification lines (IWalborn et al. 2002b) 
are not fully compatible with the helium lines traditionally 
used for classification. 

Hi) The spectroscopically determined masses of the dwarf and 
giant stars in our set of programme stars are found to be 
systematically smaller than those derived from non-rotating 
evolutionary tracks. For helium enriched dwarfs and giants, 
i.e. those having Yu e > 0.1 1, we find that all show this mass 
discrepancy. The same was f ound in an analysis of SMC 
stars using the same methods (IMokiem et al . 2006). We in- 
terpret this as evidence for efficient rotationally enhanced 
mixing leading to the surfacing of primary helium and to 
an increase of the stellar luminosity. 

iv) The bright giants and supergiants do not show any mass 
discrepancy, regardless of the surface helium abundance. 
This also is consistent with the finding for Galactic and 
SMC class I-II objects stud ied with the same methodology 
( Moki em et al . 2005, 2006). This suggests that shortly after 
birth all these stars must have rotated at less than about 30 
to 40 percent of the surface break-up velocity. 

v) A weak correlation is found between microturbulent veloc- 
ity and surface gravity. More extended atmospheres (i.e. 
lower gravity stars) require a relatively large « tur b to fit the 
lines. The reason for this relation is unclear, however, it 
does not seem to be connected to the lines being formed 
closer to the sonic point of the wind flow in low gravity 
stars. 

vi) From a comparison of modified wind momenta D mom we 
find that the wind strengths of LMC stars are weaker com- 
pared to Galactic stars, and stronger compared to SMC 
stars. Comparing the derived D mom as a functio n of lum i- 



nosity with predictions for LMC metallicities by Vink et al 



( 200 1 ) yields good agreement in the entire luminosity range 
that was investigated (5.0 < L/Lq < 6.1). 
vii) We have determined an age of, respectively, ~7. 0+1.0 Myr 
and -3.0+1 .0 Myr for the clusters LH9 and LH10. The age 
difference and relative distances are in good agreement with 
a sequential star formation scenario, in which stellar activ- 
ity in LH9 triggered the formation of LH10. 
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Appendix A: Fits and comments on individual 
objects 

The observed spectra shown in this section were corrected for 
radial velocities. If not noted differently the lines that were fit- 
ted are the hydrogen Balmer lines Ha, Hy and H<5; the Hei 
singlet line at 4387 A; the Hei triplet lines at 4026, 4471 and 
4713 A; and the Hen lines at 4200, 4541 and 4686 A. Over 
plotted are the best fit spectra, unless noted differently. For a 
discussion of the line weight ing scheme adopted in our fitting 
procedure we refer to lPaper II (especially to their Table 3.). 



N1 1-004 (Fig. |ATJ With exception of Hei ,14471 good fits 
for all lines of this OC9.7 supergiant were obtained. The under 
prediction of the He i ,1447 1 line is possibly connec ted to the 
so-called generalise d dilution effect (IVoels et all 19891 see also 
Repolust et ai]|2004 . 



N1 1-008 (Fig. IA.lt All lines are reproduced well. Note that 
given the weakness of He n /14686 and lower weight assigned 
to this line, the reproduction of its profile can be considered 
good. 



tive temperature higher by the upper error estimate we derived 
for T e ff. For even higher temperatures the He i /14471 would, in 
contradiction to the observations, disappear completely (also 
see Sect.HTTTi. 



N1 1-032 (Fig. |A2) The Hei and Hen blend at 4026 A was 
not observed for this object. In the final fit, therefore, also the 
H6 Balmer line is shown. 

N1 1-033 (Fig. IA21 This B0 giant is fast rotator. To repro- 
duce the observed line profiles a projected rotational velocity 
of 256 kms -1 was necessary. In the presented fit the strength 
of He ii ,14541 seems to be under predicted. However, given the 
relative weakness of this line this small discrepancy is negligi- 
ble. Note that the sharp change in the central part of the Ha 
profile is the result of an over subtraction of the core nebular 
feature and was not taken into account in the fitting procedure. 



N1 1 -036 (Fig.E2) The line profiles of Nl 1-036 could be re- 
produced quite accurately because of its relatively slow rota- 
tion. 



N1 1-026 (Fig. IA.lt A high effective temperature of 53 kK 
was needed to fit the spectrum of this 02 giant. Note that even 
though on its relative plot scale the very weak Hei ,14471 line 
is not fitted perfectly, the quantitative fit quality is good and 
is comparible to that of the Hen lines. Still, trying to assess 
whether the fit could be improved we ran test fits with increas- 
ing relative weight of the Hei ,14471 line. It turned out that an 
increase of a factor of five was necessary to obtain a solution 
with an improved He i ,1447 1 line fit, which is shown in Fig. IA.ll 
as a set of dotted lines. This best fit has a T e g lower by 3.7 kK 
and all other fit parameters approximately equal compared to 
the regular fit. The overal fit quality of the other diagnostic lines 
is somewhat less, though this is difficult to discern by eye be- 
cause these lines are relatively strong (one has to zoom in to the 
0.5 to 1 percent level, as we do for He i /14471). Note, however, 
that the error estimates for T e ff are relative large and within the 
lower error estimate of 3.9 kK (cf. Tab. [3]) the solutions agree. 
Also note that even for the lower T e ff solution this 02 star is 
still considerably hotter than an 03 star (see Sect. l4.lt . 



N1 1-029 (Fig. IA.lt A good fit was obtained. The slight un- 
derestimation of the core strength of the Hen ,14200 and 
He ii ,14541 lines is, given the relative weakness of these lines, 
not significant. 

N1 1-031 (Fig. |A2| ) To fit the spectrum of this 02 star a 
relatively low effecti ve temperature of 45 kK was necessary. 



Walborn et alj d2004l) derived a T eB of 55 kK for this object. 
However, their determination of this parameter was based on 
the analysis of the N iv-v ionisation adopting a fixed gravity of 
log g = 4.0. Based on the analysis of the helium ionisation we 
exclude a r e (j higher than ~47 kK. This is illustrated by the dot- 
ted lines, which correspond to a model calculated for an effec- 



N1 1-038 (Fig. IA.3t The relatively poor reproduction of 
He i /14387 in this 05 bright giant is the result of the low weight 
assigned to this line based on the spectral type of Nl 1-038. A 
good fit to the He i ,14471 could not be obtained due to its pecu- 
liar line profile shape. Possibly this triangular profile is related 
to macroturbulence. 



N1 1 -042 (Fig. IA.3t A relatively low projected rotational ve- 
locity was found for this B0 giant allowing for a nearly perfect 
fit. 



N 1 1 -045 (Fig. IA.3t A good fit was obtained for this 09 giant. 
Despite the good fit quality a problem is apparent when M s and 
M ev are compared in Tab. [2] The spectroscopic mass is found 
to be less than half its evolutionary equivalent. 

N1 1 -048 (Fig. lA~3l Based on t he fact that the He n ,14686 is 
the strongest line in its spectrum Park er et al. (1992) classified 
this star as a Vz star. The fit parameters, however, place this 
object at a considerable distance from the theoretical ZAMS in 
Fig. |5] We also find that the spectrum can not be reproduced ac- 
curately. The width of the neutral helium lines is systematically 
unpredicted, whereas the width of the He ii is over predicted. 
This and the fact that the helium abundance that is recovered 
from the spectrum is rather low (Yue = 0.06) are indicative of 
a possible binary nature. 



N1 1 -051 (Fig. |A4J| The spectrum of the 05 dwarf Nl 1-05 1 
shows very broad lines, indicating fast stellar rotation. To ob- 
tain the final fit a projected rotational velocity of 333 kms -1 
was required. 
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Fig.A.l. Comparison of the observed line profiles of Nl 1-004, -008, -026 and -028 with best fitting synthetic line profiles 
obtained using the automated fitting method (grey lines). Wavelengths are given on the horizontal axis in A. The vertical axis 
gives the normalised flux. Note that this axis is scaled differently for each line. The dotted line profiles for Nl 1-026 correspond 
to the best fit obtained adopting a five times larger relative weight for the He i A4471 line. This best fit has a T e ff lower by 3.7 kK 
and all other fit parameters approximately equal compared to the fit for the nominal He i A4471 weight. 



N1 1-058 (Fig. IA.41 > A good fit was obtained for all lines of temperature we find for these two objects are also compara- 

this 05.5 dwarf. The small under prediction of the Hei4387 ble, the helium spectrum suggest a spectral type of 03 for both 

line strength is the result of the relatively low weight assigned objects, 
to this line for the spectral type of this star. 



N1 1 -060 (Fig. |A4l! Note that the He i ,14471 line of this 03 
dwarf has a strength that is comparable to the strength of this 
line in the spectrum of the 02 giant Nl 1-031. As the effective 



N 1 1 -061 (Fig. \A.4\ All lines of this 09 dwarf are reproduced 
correctly. The mass loss could be reliably determined using the 
automated method at a rate of 2.1 x 10~ 7 Mq. Note that the 
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Fig. A.2. Same as Fig.|AT] however, for Nl 1-031, -032, -033 and 
a 2.2 kK increase in r e ff . 



036. Shown as dotted profiles for Nl 1-031 is the effect of and 



value of 1.8 obtained for the wind acceleration parameter ft is N1 1-066 (Fig. IA.5b Apart from a slight under prediction of the 



rather high for a dwarf type object. 



width of He n /14686 all lines of this 07 dwarf could be fitted 
with good accuracy. Like the previously discussed object, again 
a large mass discrepancy is found with a ratio of M ev to M s of 
1.8. 



N1 1-065 (Fig. IA.5b A relatively high helium abundance of 
0.17 was recovered from the spectrum. The good fit quality 
obtained for all lines, however, does not indicate an overesti- 
mation of this parameter. A strong mass discrepancy is found 
with M s being smaller than M ev by approximately a factor of N1 1-072 (Fig 
two. 



N1 1 -068 (Fig. IA.5b The final fit reproduces the observed line 
profiles very accurately. No further comments are required. 



The star has a very sharp lined spectrum. 
The best fit required a v T sin / of 14 km s -1 . Note that given the 
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Fig. A.3. Same as Fig.|AT] however, for Nl 1-038, -042, -045 and -048. 



lower error estimate, given in Tab. [3] and the spectral resolution hire needed to fit the spectrum is, however, higher by more than 
of the data this value can also be interpreted as an upper limit. 2 kK. This is the result of a log g that is higher by ~0.2 dex. 



N1 1 -087 (Fig. lAT6l > To fit the spectrum of Nl 1-087 a v r sin i 
of 276 kms 1 was needed. The final fit shows a good repro- 
duction of the line profiles with a slight under prediction of the 
He ii /14541 line. Given the signal-to-noise ratio of the spectrum 
and the relative strength of this line, we do not believe this to 
be significant. 



N1 1 -123 (Fig. IA.6l ) Like the previously discussed object Nl 1- 
123 is also classified as an 09.5 dwarf. The effective tempera- 



Bl 237 (Fig. IA.7I ) A good fit was obtained for this 02 dwarf 
star. The He 1 1447 1 line is hardly visible. Fortunately, the hy- 
drogen and Hen lines provide the automated fitting method 
with enough information to determine an effective temperature. 
The error estimates on this parameter, however, are consider- 
able (see Tab. 0. The effect of lowering T e ff to the lower error 
estimate of 49 kK is illustrated by the dotted profiles. 



BI 237 was recently also analysed bv lMassev et al. (2005). 
Compared to this study we find large differences for T e g 
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Fig. A.4. Same as Fig.|AJ] however, for Nl 1-051, -058, -060 and -061. 



(+5 kK), logg (+0.2 dex) and M (-0.4 dex). Our solution 
for higher effective temperature and surface gravity can be ex- 
plained by the fact that the ionisation structure of_the atmo- 
sphere is set by both these parameters (also see IPaper II) . The 
reduced mass loss rate we find is the result of the increased 
value for f3 (+0.4) obtained by the automated method. 



Bl 253 (Fig. IA.7b The star has as an identical spectral type 
as the previously discussed object. The photospheric parame- 
ters determined from the spectrum are also very similar to the 
parameters of BI 237. Note that in contrast to this agreement 
BI 253 is found to have a much denser wind. Its mass loss rate 



is larger by a factor 2.5 and reflects the fact that compared to 
BI 237 the wind lines Ha and He n /14686 are more filled in. 



In their analysis iMassey et alJ (120051) could only determine a 
lower estimate for T e g of 48 kK, which is consistent with our 
lower error estimate. The other parameters determined by these 
authors are in fair agreement with our analysis. An e xception 



estimate 



to this is the surface gravity for which IMassey et al. 
a value lower by 0.3 dex, which is the result of the lower T e g 
used for their fit. 



M. R. Mokiem et al.: Wind properties and evolution of hot massive stars in the LMC 



25 




N1 1-065 06.5 V((f)) 

Ha Hel + Hell 4026 



1.00 S£T ' "u^*f_ 



0.80 - 




4330 4340 4350 6547 6562 6577 4021 4026 4031 

Hel4387 Hel4471 Hel4713 

1.00 



0.94 - l ' 30 ' 



1.00 



0.90 - 




4383 4388 4393 4468 4471 4474 

Hell 4200 Hell 4541 



4711 4713 4715 
Hell 4686 




1.00 



4195 4200 4205 




4537 4541 4545 



4682 4686 4690 



N1 1-066 07V((f)) 

Ho 



Hel + Hell 4026 




4330 4340 4350 6547 6562 6577 4021 4026 4031 

Hel 4387 Hel 4471 Hel 4713 

1.00 



0.96 - \ P ~ 

o.92 - \r - 




4385 4388 4391 4468 4471 4474 

Hell 4200 Hell 4541 



4711 4713 4715 
Hell 4686 




1.00 




4195 4200 4205 



4537 4541 4545 




4682 4686 4690 



1.00 



0.80 - 



1.00 
0.90 - 
0.80 - 

0.60 ■ '■ ] 0.70 - 

4330 4340 4350 6547 6562 6577 

Hel 4387 Hel 4471 




N1 1-068 07V((f)) 

Ha Hel + Hell 4026 

1.00 





4021 4026 4031 
Hel 4713 




4385 4388 4391 4468 4471 4474 



4711 4713 4715 
Hell 4686 





41 95 4200 4205 



4537 4541 4545 



4682 4686 4690 



N1 1-072 B0.2III 

Ha 



Hel + Hell 4026 




4386 4388 4390 
Hell 4200 



4468 4471 4474 4712 4713 4714 





4192 4200 4208 



4536 4541 4546 



4684 4686 4688 



Fig. A.5. Same as Fig.|AT] however, for Nl 1-065, -066, -068 and -072. 



Sk -66 18 (Fig. IA.7I ) With exception of the Hen lines at per error estimate. The effect of this increase of 0.12 dex on the 
4541 A and 4686 A we obtained a near perfect fit. The helium line profiles is shown using dotted lines, 
abundance Tne = 0.14 indicates an enriched atmosphere. 



Sk -6 6 1 00 (Fig. IA.7t was previously analysed bv lPuls et al 



Compared to this study we find similar parameters with 
exception of a lower T e s which is the result of the inclusion of 
blanketing and an increased helium abundance of Fee = 0.19. 
The final fit shows that with this large value all lines could 
be reproduced quite accurately. The slight mismatch of the 
He n ,14686 line can be improved by increasing M to its up- 



Sk -67 166 (Fig. IA.8b This supergiant has both Ha and 
He n /14686 in emission indicating the presence of a dense stel- 
lar outflow. To correctly reproduce the line profiles a very high 
mass loss rate of 9.3 x 10~ 6 Mgyr -1 was required. This value is 
in good agreement with the findings of ICrowther et al.1 (|2002) 
who s tudied this objects usin g the model atmosphere code cm- 
fgen dffillier& Miller! Il998l) in the optical, UV and far-UV. 
Compared to these authors we also find that the other fit pa- 
rameters are in good agreement. An exception to this is the he- 
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Fig. A.6. Same as Fig.|AJ] however, for Nl 1-087 and -123. 
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Hum abundance. ICrowther et al.l d2002l) adopted a fixed value 
of ^He = 0.2 whereas our automated fitting method was able to 
self consistently determine a value of Yn e = 0.28. 



Sk -70 69 (Fig. |AJ) To fit the spectrum of this 05 dwarf 
Sk -70 69 a helium abundance of Yu e = 0.17 was required. 
Note that, similar to Nl 1-065, which is also strongly enriched 
with helium, the spectroscopic mass is found to be much 
smaller than the evolutionary mass. 
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Fig. A.7. Same as Fig. |AJ] however, for BI 237, BI 253, Sk -66 18 and Sk -66 100. The dotted profiles for BI 237 correspond 
to a model calculated with a T e s reduced by 3.8 kK compared to the best fit value. For Sk -66 100 the dotted profiles show the 
effect of an increase in M by 0.12. 
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Fig. A.8. Same as Fig.|Aj] however, for Sk -67 166 and Sk -70 69. 
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